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ABSTRACT 

The  non-specular  reflection  of  plane  electromagnetic  waves  of  arbi- 
trary polarization  by  certain  perfectly  conducting  surfaces  composed  of 
either  semi cylindrical  or  hemispherical  bosses  on  an  infinite  plane  is 
analyzed.  Solutions  in  terms  of  eigenfunctions  for  the  problem  of  the  single 
boss  on  an  infinite  plane  and  a  plane  wave  at  an  arbitrary  angle  of  incidence 
are  derived  and  extended,  subject  to  the  6ingle  scattering  hypothesis, to  ob- 
tain the  far  field  solutions  for  certain  small  finite  patterned  distributions 
and  both  small  finite  and  infinite  uniform  random  distributions  of  bosses 
small  compared  to  the  wavelength.  The  results  for  the  various  cases  are 
then  compared  in  the  plane  of  incidence  and  similarities  between  the  analogous 
expressions  for  the  distributions  of  semicylinders  and  hemispheres  ere  noted. 
Expressions  are  obtained  for  the  ratios  of  the  reflected  intensities  and 
radial  energy  flux  polarized  parallel  and  perpendicular  to  the  plane  of  inci- 
dence's well  as  for  the  total  intensity  and  radial  energy  flux  for  the 
case  where  the  incident  wave  is  unpolarized.  It  is  found  that  for  certain 
values  of  the  parameters  the  reflected  radiation  may  consist  only  of  either 
the  specular  or  the  scattered  contributions,  while  for  other  values  of  the 
parameters  one  of  the  scattered  contributions,  either  parallel  or  perpendi- 
cular component  may  vanish.  The  results  also  indicate  the  occurrence  of  an 
extremum  in  the  reflected  radiation  in  the  vicinity  of  the  specular  angle  of 
reflection  which  for  certain  ranges  of  the  parameters  for  the  small  finite 
distributions  may  be  a  minimum  rather  than  a  maximum.  For  these  cases  there 
is  also  some  critical  angle  of  incidence  (not  necessarily  Tf/2   or  grazing 
incidence)  for  which  the  reflection  at  the  6pecular  angle  is  completely 
specular. 

The  analogous  distributions  of  cylinders  and  spheres  are  also  con- 
sidered. 
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NOTITIOH 

E,  H,  -  The  electric  and  magnetic  field  intensities  which  comprise  the  electro- 
magnetic field . 

i  -  Superscript  or  subscript  character! ?ing  the  Incident  wave, 

r  -  Superscript  or  subscript  character! 7ing  the  total  reflected  wave. 

p  -  Superscript  or  subscript  characterizing  the  specularly  reflected  rave  components, 

c,  8  -  Superscript  or  subscript  characterizing  the  scattered  components  of  the  total 
reflected  wave  for  the  (semi) cylindrical  and  (hemi) spherical  cases  re- 
spectively. 

A,  B,-  Arbitrary  constants  (determined  from  proscribed  initial  conditions)  that 
specify  the  polarization  parallel  and  perpendicular  to  the  plane  of 
incidence. 

II',  II"  -  Vector  Hertz  potentials  of  the  electric  and  magnetic  type  respectively. 

y^yw  _  scalar  components  of  the  Hertz  potentials, 

K  =  Kn  =  (27T/71   )n  =  (ou/cjfi"  -  Propagation  vector  (the  symbol^  denotes  the 

unit  vector) . 

o(,Q     ,  -  Angles  of  incidence, 

i,i,k  -  Unit  rectangular  coordinate  vectors, 

R  =  r(cos  9  i  +  sin  8  j)+  zk  -  Position  vector  in  cylindrical  coordinates  . 

A1,  B*  -  Constants  introduced  for  brevity. 

a  -  Radius  of  the  boss. 

£  n  -  Neumann's  factor:  €  n  -  1  if  n-0,  £n  =  2ifn>0 

J   -  Bessel  functions  , 
n 

B_  -   Hankel  functions  of  the  first  kind, 

Zn  -  Derivatives  of  the  cylindrical  functions  with  respect  to  their  arguments. 

A  ,  A  -  Cylindrical  scattering  coefficients, 

o  -  Left  superscript  denoting  the  approximate  solution  (Kr  »  1,  Ka < 1)  for  the 
single  boss  at  the  origin  of  the  reference  coordinates. 

Q_  (or  K) ,  c i  -  Abbreviations  for  the  amplitude  and  phase  of  the  scattered  com- 
ponents. 

F  -  Either  E  or  H 
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a  -  Summation  index  ranging  from  -N  to  +  N  , 
?!)  -   Total  number  of  bosses  in  e  distribution, 

/  -  Fraunhofer  phase  factor  , 

j  -  Subscript  labeling  the  superposed  distributions. 
/?? ,  -  Factor  defined  so  that  'V .  -   0  for  the  particular  value  associated 
with  the  bosses  of  largest  radii  and  *y .  -   1  for  all  other  values. 

M  -  Density  function  eoual  to  the  number  of  bosses  per  square  centimeter  • 

d  -  Linear  extent  of  the  small  finite  distributions  • 

C,  D  -  Constant  coordinates  of  a  field  point  for  the  infinite  distributions, 

V  -  Volumetric  departure  from  the  plane  equal  to  the  total  volume  of  bosses 
per  square  centimeter  of  distribution. 

r  =  r  (sin  9  cost^  i  t  sin  6  eir.tf    J  +  cos  9  k)-  Position  vector  in  spheri- 
cal coordinates. 


-  Spherical  Bessel  functions  -  Morse's  notation  , 


D  ,  etc.  -  Partial  derivative  with  respect  to  r,  etc. 
r 

nnm'  Gnm»  ^  Pnm»  Pnm  *  Symbols  introduced  for  brevity* 

P  (cos  9)  -  Derivatives  of  the  Associated  Legendre  functions  with  respect  to  9, 

Xn5m  ,  2£  n^-m  -  Sums  over  n  and  m  such  that  n  and  m  have  odd  or  even  parity 
respectively^ 

a  ,  a  -  Spherical  scattering  coefficients  . 
n   n 

y   -  Angle  made  by  a  line  of  boeses  with  the  x-axis  in  the  xy-plene:  also 
used  to  label  a  particular  line  of  bosses. 


O     ,u      -   Polar  coordinates  in  xy-plane 

Angle  made  by  the  axes  of  a  recta 

the  xy-axes:  also  used  to  label  a  particular  rectangle. 


S  -  Angle  made  by  the  axes  of  a  rectangular  distribution  of  spheres  with 
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u  -  "Reflection  coefficient"  for  the  infinite  distributions, 

Y     -  Phase  factor  of  u  . 

A  -  function  -  A  function  such  as  sin  P  /P  or  sin  OOP  /sinP  which  has  a 
single  principle  maximum  or  a  series  of  such  maxima  and 
which  reduces  to*?7  at  the  specular  angle  of  reflection  , 

J  s  Re(E«E  )  -  Quantity  proportional  to  the  intensity 


S  -  Padial  component  of  the  complex  Poynting  vector, 

<r  =  p„/sx 
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1.  Introduction 

The  scattered  reflection  of  plane  scalar  waves  from  certain  sur- 
faces consisting  of  distributions  of  serticylindrical  or  hemispherical 
bo-ses  (protuberances)  on  an  infinite  plane  was  analyzed  in  a  previous 
report  .  The  present  report  considers  the  analogous  problems  for  the  re- 
flection of  electromagnetic  waves  from  metallic  surfaces  of  infinite  conduc- 
tivity. 

The  procedure  we  will  follow  will  be  to  obtain  solutions  in  terms 
of  the  appropriate  eigenf unctions  for  the  reflection  problem  of  a  plane  elec- 
tromagnetic wave  of  arbitrary  polarization  incident  at  some  arbitrary  angle 
on  either  a  semicylindrical  or  heni spherical  boss  on  an  infinite  plane.  The 
solution  for  the  single  boss  will  then  be  extended, subject  to  the  "single 
scattering"  hypothesis, to  obtain  approximate  far  field  solutions  for  certain 
small  (as  compared  to  the  distance  of  observation)  finite  patterned  distri- 
butions and  both  small  finite  and  infinite  uniform  random  distributions 
of  bosses  of  radii  small  compared  to  the  wavelength.  A  discussion  of  the  de- 
rived results  will  be  reserved  for  the  last  section  where  the  analogous  ex- 
pressions for  the  cylindrical  and  spherical  cases  will  be  compered  in  the 
plane  of  incidence. 

Except  for  the  basic  solutions  to  the  problem  of  the  single  boss,  the 
mathematical  scope  and  content  of  this  report  do  not  differ  appreciably  from  the 
previous  treatment  of  the  scalar  cases.  For  this  reason,  while  the  derivation 
of  the  solutions  to  the  problems  of  the  single  boss  will  be  given  in  detail, 
only  the  significant  differences  between  the  present  and  previous  procedures 
for  deriving  the  results  for  the  distributions  will  be  included—the  manipu- 
lations required  in  going  from  the  sinple  boss  to  it-«?  distributions  being  very 

1.  V.  Twersky,  On  the  Scattered  Reflection  of  Scalar  Waves  from  Absorbent 
Surfaces,  New  York  University,  Mathematics  Research  Group,  Research 
Report  No.  EM-??,  Aupust,  1950. 


similar  to  those  employed  oreviously  and  leading  eventually  to  essentially 
the  same  sums  and  integrals  which  were  there  evaluated.   (Page  references  to 
the  previous  report  will  be  found  in  the  text  to  facilitate  the  discussion). 
The  eigenfunction  formulation  employed  for  the  problem  of  the  single 
boss  restricts  us,  at  least  for  practical  purposes,  to  the  range  where  the 
radii  of  the  bosses  are  small  compared  to  the  incident  wavelength  (or  more 
strictly  speaking  where  Ka  =  2  IT  a/>.  <  1).  Although  the  functions  employed 
are  tabulated  and  in  principal  the  initially  derived  expressions  are  valid 
for  all  values  of  Ka,  it  is  only  for  Ka<  1  that  the  serieB  are  sufficiently 
convergent  to  be  approximated  by  their  first  few  terms.  The  expressions  for 
the  distributions  derived  subject  to  this  limitation  sre  therefore  only  ap- 
propriate for  the  diffraction  range  of  phenomena  (or  the  "Rayleigh  scattering" 
region)  as  might  occur  with  surfaces  whose  irregularities  are  small  compered 
to  the  wavelength.  In  this  region  the  distributions  of  bosses  on  a  plane 
present  features  characteristic  of  either  striated  or  rough  surfaces  and 

should  be  useful  in  investigating  some  of  the  general  aspects  of  non-specular 

2 

reflection. 

The  restriction  of  Ka<  1  of  necessity  precludes  the  discussion  of  a 
wide  variety  of  phenomena.  In  particular  it  excludes  the  optical  range 
where  there  is  reason  for  believing  that  reflection  is  governed  by  phenomena 

arising  because  the  usual  surface  irregularities  are  large  compered  to  wave- 

2.  As  was  mentioned  in  the  previous  report,  the  solution  for  the  simple  semi- 
cylindrical  boss  had  been  obtained  by  Rayleigh,  reference  6.  It  is  in- 
teresting to  note  that  another  approach  to  such  problems  as  reflection 
from  a  rough  surface  introduced  by  Rayleigh, "Theory  of  Sound",  Section 
2728,  Dover,  N.T.,  (1945) — who  considered  the  problem  of  a  plane  wave 
incident  on  a  corrugated  surface  whose  cross-section  could  be  repre- 
sented by  a  Fourier  series — is  being  developed  by  C.  T.  Tai,  "Reflection 
and  Refraction  of  a  Plane  Electromagnetic  Wave  at  a  Periodical  Surface", 
Harvard  Tech.  Report  No.  28,  Cruft  Lab.,  (Jan.  194-8)  and  S.  0.  Rice,  "Wave 
Propagation  Over  Rough  Surfaces",  paper  delivered  nt  the  Symposium  on 
the  Theory  of  Electromagnetic  Waves,  New  York  University,  June,  1950. 
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length — essentially  the  region  for  which  the  empirical  Lambert '3  Law  of 
diffuse  reflection  is  considered  valid.  A  certain  amount  of  theoretical 
work  has  been  done  in  this  range,  particularly  in  an  effort  to  derive 
Lambert's  Law  which  is  believed  to  be  a  statistical  effect  of  specular  re- 
flections occurring  at  the  irregularities  themselves.-5  This  phase 
of  the  problem  will  perhaps  be  investigated  by  the  techniques  of  geometrical 
optics  in  a  future  report. 

The  region  Ka<l,  however,  is  one  in  which  some  interesting  phenomena 
are  predicted.  It  is  here  found,  for  example,  that  for  certain  ranges  of  the 
parameters  the  results  are  at  variance  with  what  might  be  surmised  from  an 
intuitive  basis  and  yield  a  minimum  for  the  reflected  radiation  observed  in 
the  vicinity  of  the  specular  angle  of  reflection.   They  also  predict  the  oc- 
currence of  certain  critical  angles  of  incidence  for  which  the  reflection  at 
the  specular  angle  is  completely  specular,  and  certain  critical  angles  for 
which  marked  changes  occur  in  the  ratio  of  the  reflected  components  nolarized 
parallel  and  perpendicular  to  the  plane  of  incidence.  It  would  be  of  interest 
to  subject  this  region  (Ka<  1)  to  a  critical  experimental  investigation  to 
determine  whether  or  not  these  predictions  based  on  relatively  simple  models 
can  be  verified  for  actual  surfaces. 

2 .  Formulation  of  the  Single  Boss  Problem  and  Some  Preliminary  Considerations 

Consider  a  plane  electromagnetic  wave  of  arbitrary  polarization  in- 
cident at  an  arbitrary  angle  on  a  perfectly  conducting  surface  consisting  of 
either  a  semi cylindrical  or  a  hemispherical  boss  on  an  infinite  plane.  We 


3.  W.  W.  Barkas  and  R.  F.  S.  Hearman,  Phys.  Soc.  Proc,  51  274  (1939);  W.  A. 

Rense,  J.  Opt.  Soc.  Amer.,  40  55  (1950);  L.  Oraestine  and  P.  Van  der 

Berg,  Physica,  4  1181  (1937)7 
4..  Such  minima  at  the  specular  angle  of  reflection  were  actually  observed  for 

the  acoustical  case  with  sound  incident  on  surfaces  whose  irregularities 

were  of  the  same  order  of  magnitude  as  the  wavelength.  See  E.  Skudrzyk, 

Akustische  Zeits.,  4  189  (1939). 
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seek:  solutions  of  Maxwell' 6  equations  that  satisfy  the  boundary  condition  that 

the  tangential  component  of  the  electric  field  Intensity  vanishes  at  the  surface, 

and  are  subject  to  the  usual  Sommerfeld  radiation  condition  that  their  scattered 

components  reduce  to  outgoing  waves  at  infinity.  The  incident  polarization 

vector  or  electric  field  intensity,  B  ,  is  resolved  into  components  parallel  and 

perpendicular  to  the  plane  of  incidence  (AE  ,  BE  ;  the  constants  so  introduced 

to  be  determined  from  prescribed  initial  conditions)  and  we  seek  to  determine  the 

reflected  field  specified  by  E_r  and  Hr  such  that  B  B   B  ♦  Er  and  H  *  H  ♦  Hr 

satisfy  Maxwell's  equations,  the  boundary  condition  B   I    -    =0,  and  the 

radiation  condition. 

The  problem  of  scattered  reflection  by  the  boss  on  the  plane  is  soluble  by 

an  image  method  whereby  the  original  problem  is  replaced  by  that  of  the  scattering 

by  either  a  cylinder  or  sphere  of  two  incident  waves — the  original  wave  specified 

by  E  and  its  image  in  the  removed  plane  specified  by  E;  (where  Er  is  determined 

from  the  boundary  condition  (Et   +  S?  )  ,    =  0) — and  the  determination  of 

tan    tan  (plane 

the  scattered  field  EC*8.  The  configuration  of  Fig.  la  (or  Pig.  2a)  Is  replaced 
by  that  of  Pig.  lb  (or  Pig.  2b)  and  we  proceed  to  solve  the  image  problem,  since  it  can 
be  shown  that  the  solutions  in  regions  bl  are  identical  with  those  in  regions  al,  the 
required  solutions,  satisfying  all  the  required  conditions:  I.e.,  that,  for  example, 
the  electric  field  intensity  in  la,  E  »  E  *  Er,  is  identical  with  that  in  lb  for 
|8|^  7T/2,  E  *  S  *  E?  ♦  EC,  so  that  the  determination  of  B  (the  component  scattered 
by  the  cylinder  obtained  by  imposing  the  boundary  condition  BtaIJCyiinder  =  °»  and 
so  chosen  so  that  it  also  satisfies  the  required  boundary  condition  on  the  plana) 
will  yield  the  scattered  reflected  component  for  the  boss  on  the  plane,  B  =  B^  ♦  E  . 

Wa  will  follow  the  above  procedure,  but  it  will  be  found  more  convenient 
to  formulate  the  problems  in  terms  of  suitable  Herti  vector  potentials  from 
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which  the  field  may  be  resiily  derived.  It  can  be  shown  that  the  elec- 

_iwt 
tromagnetic  field  in  free  space  for  a  harmonic  time  factor  e     .  E,  H, 

may  be  resolved  into  two  partial  fields  E  ,  H  and  E  ,  H  ,  derivable  from 
the  Hertz  vector  potentials  II',  11%  where  the  primed  variables  are  of  the 
"electric  type"  (associated  with  a  distribution  of  electric  dipoles)  and  the 
double  primed  variables  of  the  "magnetic  type"  (associated  with  a  distribution 
of  magnetic  dipoles).  In  mks  units  the  relationship  between  the  above  quan- 
tities may  be  expressed  as 

(1)  E«   =     VxVxII'   =   (K2      +W-)nt  f       hi   =     £7x11  •   =-iKyV  xW  f 

(2)  E"  =     -ptf  xll"  =  iyujV  xll"         }       H"  =VxVxII"  =   (K2+W  .)H"    , 
where  II'   and  II"  are   solutions  of  the  vector  wave  equation, 

(3)  VxVxII  _  WTI+    /^£    II     =     (VxVx-W*     -K2)II  =  0  , 
JV  =   (27T/A)2  =   (W/c)2  =C^u.2t  ,y  =  K/ffW      =   (€/^)lj 

where  II  in  (3)  stands  for  either  II*  or  II".  (For  brevity  and  where  no  con- 
fusion will  arise  ve  will  also  in  the  future  use  II  to  denote  either  II'  or 
II"  and  similarly  with  their  associated  constants). 

Tor  the  problems  with  which  we  will  be  concerned,  the  scalar  compon- 
ents of  II  will  consist  of  sets  of  either  the  cylindrical  or  spherical  eigen- 
functions  with  arbitrary  constants  to  be  determined  from  the  boundary  and 
initial  conditions.  Tor  the  cylindrical  case  it  will  be  convenient  to 
choose  the  Hertz  vectors  of  the  form  II  ■  j  ic,  parallel  to  the  axis  of  the 
boss  along  z,  while  for  the  spherical  ease  we  will  choose,  II  ■  V"r,  in  the 
radial  direction.  Tor  these  choices  of  XI  it  will  then  be  recognized  that 
the  primed  and  double  primed  quantities  in  equations  1-3  are  tne  perhaps  more 
familiar  TM  (transverse  magnetic)  and  13  (transverse  electric)  modes  re- 

5,  See  for  example,  J,  A.  Stratton,  Electromagnetic  Theory,  pp.  28-32,  McGraw- 
Hill,  H.  Y.,  19U1. 
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spectivelyj  the  term  transverse  meaning  perpendicular  to  the  vectorial 
direction  of  II  and  not  necessarily  to  its  direction  of  propagation,  i.e., 
thus  II'  indicates  that  its  associated  magnetic  field  has  no  component  parallel 
to  k  for  the  cylindrical  case  or  no  component  parallel  to  r  for  the  spherical 
case.  In  the  following  sections  which  are  concerned  with  the  mathematical 
analysis  of  the  problems  it  will  be  seen  that  equations  1-3  simplify  consider- 
ably because  of  this  choice  of  the  II. 

3.  The  Semicylindrical  Boss. 

A.  The  Single  Boss  . 

Consider  a  plane  electromagnetic  wave  characterized  by  jir  incident 
on  a  semicylindrical  boss  of  radius-a  with  axis  parallel  to  the  z-axis  on  an 
infinite  y7 -plane  as  in  Fig.  la.  We  seek  a  solution  of  the  problem (as  formu- 
lated in  2  )  subject  to  the  boundary  condition  that  the  tangential  component 
of  E  vanishes  on  the  plane  ani  boss  or  that 

U)        (^--1^  =  0)^^  ,(^  =  ^  =  0)^  . 

The  image  technique  discussed  in  the  preceding  section  will  be  employed  to 
solve  this  problem. 

In  general,  the  incident  field  can  be  initially  decomposed  into  com- 
ponents parallel  and  perpendicular  to  the  plane  of  incidence (as  determined  by 
the  unit  vectors 'n.  and  i),  expressed  in  terms  of  the  reference  cartesian  co- 
ordinates and  then  transformed  to  the  cylindrical  coordinates  of  the  boss,  as 

follows  sequentially: 

6.  This  technique  was  employed  previously  by  Rayleigh,  Phil.  Mag.,  1£  350 
(1907)  to  solve  the  above  problem  for  the  case  of  incidence  perpendicular 
to  the  axis  of  the  boss. 


Figure  la 


Figure   lb 


/  plane   of 
incidence 


a. 


b. 


II  (region   of  X  negative)  I  (region   of  X  positive) 

A  plane  wave  incident  on  a  semi cylindrical  boss  (whose  axis  is  parallel  to  the 
z-axis)  which  lies  on  an  infinite  yz-plane.  The  plane  of  incidence  is  the  plane 
of  ILi  and  the  z-axis.  The  unit  vector  ni  nakes  an  angle  £  with  the  z-axis  while 
Its  trace  on  the  xy-plane  makes  an  angle  <*  with  the  x-axis. 

Two  plane  waves,  n\  and  its  image  in  the  y*-plane,  lip,  simltaneously  incident 
on  a  circular  cylinder  of  radius-a  whose  axis  is  parallel  to  the  z-axis. 
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/\  S}\         S>  /N  /?• 


i  1  i  C         (n4x  l)x  n4  n.xi   )    iK, -P 


s j_\  a)  (1-sin  o    cos^  )i-sin*a   cose*  sina.  j+sino  cosfl  cos<*  k 

(l-sin^cos2*  )*  J   1  p  p  r       r  J 


-»•  B  £  cosq  j  +■  sina  sin**  k  J  > 


(1-sin  A   cos**   ) 


W    (Afcos  0  -  sin2fl  cos<x  cos   (8  -<^)J+  B  cos  a  sin  9)   r 


(5)  f  (A[-sin  9  +  sin2fl  cosotsin  (9  -crt)  +  B  cos o  cos  9)   9 

+   sin  a  (A  cosd  cos*  +  B  sin*  )   k     / 


and 


4  -x         %'S 

Hx   = 


=   /*?    n  X  F1  =/?M    -B  nix  l)yn,+  An.x  if r r — -r 

A>w  /      i    -        /<  *         *       Ml-sin2*  cos2*)? 


=    — *L-£ S     -B(l-sin2fl  cosot  jit  (Bsinfe    cos<*  sin**  A  cosr?  )  ^ 

(l-sin2o   cos2*  )H 

+  (-B  sino  cosfl  cos*  +  A  sinft  sino<  )   k  f 

iK  «P 

=  ■  *y.  e-  ~ 4    (-bFcos  9  -  sin  a   coso«  cos  (9  -d   )]  +  A  cosg  sin  9)  "r 

(l-sin2|?  cos2*  )H  l  ^  r 

f  (-pi-sin  9  +  sin2g  coso(  sin(9  -<*))-t-  A  cosfi  cos  9)   9 
•+sin  ft  (-B  cosfl  cos<A-t-  A  sin  *)   k  £  . 
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In  these  expressions  the  position  vector  of  a  field  point  P  is  R  = 
r  (cos  ©  i  ■*■  sin  9  J)  *  zk  =  r  ♦  zk;    the  incident  propagation  vector  is  K     =  Kn     ■ 
K(-sinfi.  cos*  i-sina  sinw'J  ♦  cosa  k);   and  the  transformations,   i  =  cos  6  r  -  sin  9  9. 
J  =  sin  9  r  ♦  cos  ©  9,  k  =  k,  have  been  employed.      It  can  be  seen  that  H1  may  he 
obtained  from  fi    by  replacing  A  by  -B  and  B  by  A  and  multiplying  the  resulting  ex- 
pression by/*?  . 

The  image  wave  in  the  yz-plane  (characterized  by  JSP  and  Mp)  is  the   specularly 
reflected  wave  that  would  be  obtained  were  the  boss  absent  and  is  determined  by   im- 
posing the  first  of  conditions   (U)  on  2     *  EP.     Thus  on  replacing  K,   by   the  image 
propagation  vector  K     ■  Kn"    ■  K(sing  coso(/J-8infi  sintf'j  ♦  costflcj  we  obtain  aP  by 
changing  the  signs   of  the  j  and  k  components   of  B  ,  and  Hp  by  changing  the  sign  of 
the  i  component  of  H  .     Hence  we  obtain 
i£  .a 

&  ■ — 5 = ttz)   A(l-sin2g  C082<X  )i  ♦  (A  sin2,s  cosrfeln     -Beos  6)t 

(1-sin2^  cos2*  )1/2  C  "  I  *'  J 

-(A  sinp  cosfl  cos**  B  sin©  sin*  )  Ic} 


(6) 


IK   .R 


= - — = 5 rrr  )    (Ajcos  9  -  sin  a  cos*  cos  (9  ♦  o*  jl-  B  cos/?  sin  9) "? 

(l-6in2a  cos2*  )1/2  (        L  P  e 

♦  (A  [-sin  9  *  sin2^  cos*  sin  (6  ♦  <*  j]  -B  ccsft  cos  0)  e" 

-  sing  (A  cos^  cos<^     ♦  B  sin*  )  1c? 


and 


iK  «R 


Hp  <=/*  n  x  &=     1?    e  o    "    o TTZ  S  B(l-sin2ff  cos2*  )i  ♦  (B  sin2*  cob*  sin*  ♦  A  costf  ) 

~       /     V    ~        Ulsin2^  coa2*  )1/2  L  P  P 

♦  (-B  sinfi  cosff  coso<+  A  sin/jsintf)  k£ 

iK   .R 

=  2Z 2 - -j-j  (bFcos  9  -  sin2p  cos  «  cos  (9  +  ot^  +  A  cos^sin  9)  r 

(l-sin  ft   cos*    )  '     ^ 

♦  (b[-  sin  9  ♦  sin2A  cos*  sin  (6  ♦  <*  )J    ♦  A  cos^  cos  e)  6 

♦  sine  (-B  coso  cos«<  ■►  A  sin*  )  k> 
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c      c 

To  proceed  further  and  discover  the  fields  E  and  H  scattered  by 

the  cylinder  subject  to  the  boundary  conditions  at  r  =  a,  we  no**  find  it  con- 
venient to  introduce  two  Hertz  potentials  of  the  form  II  =^k  =  f(r,  G)  eihz  k. 
Then,  since  (Vx7x  -VV-rf'k  -  -^2^,    equation  (3)  reduces  to  (72-rK2)'V'=  0, 
the  scalar  *ave  equation  for  a  harmonic  time  dependence,  so  that  the  7  are 
the  usual  cylindrical  eigenf  unctions.   We  therefore  obtain  from  (l)  and  (2) 
that  the  field  can  be  expressed  as  the  sum  of  the  two  partial  fields* 

E*   -    (K2+7<7.)CV"k>  K2^Vk+ihVr=Wih|7   VS^Is    tk(K*-h2)]  Y'^ 
(7) 

H'   =  -  iK^Vx    (t^)   =  iK^X7f*=  iK^  -r  i  |^    +3     \^\f%         , 

and 

E"  =  ipwvx   CV"k)  =  iK^I^-le   -  ^  |t?T  " 

(8) 

H"  =   (K2-r7^7  -)(y"k)  =  i  Sin  i-   +3  iil  i-     ^(K2-h2)^"        ; 

where  h  =  K  confl  and  where 

E'  -  E'1  +  E'p  +  E'C      ,       E"  -  E"1  +  E"P  +  E"c    j 

(9)    H'  =  H'1  +  H'P  +  H'c      ,       H"  =  Hpi  +  H"P  -»■  HnC 

will  now  be  determined. 

It  should  be  noted  from  (7)  and  (8)  that  X  "I  ■  (I2-  h*)^'  and 
H  »  H  ■  (K  -  h  )r".  The  following  procedure  hinges  on  these  two  expressions 
since  the  s  components  of  the  incident  and  image  waves  yield  the  corresponding 
V'  *s  with  which  we  determine  the  scattered  V'»  and  hence  the  scattered  field. 


7.  A  discussion  of  the  functions  end  expansions  used  in  this  section  will  be 
found  in  Stretton  (referenced),  Ch.  VI.   Eauations  (7)  and  (8)  will  be 
found  on  pp.  360,  361. 
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The  Y" ' .   and  *V*     are  readily  determined  from  the  z  comocnents  of   (5) 
i  p 

and   (6).     Thus   (since  K2  -  h2  =  K^    (l-cos2(3    )  =  K2  sin2  (3    )   we  obtain 
V"   =  E*/K2sin2(3    =  B*   e1^* 

V '   =  EP/K2sin2«     =-B'   e^P*£ 
p         z  ^ 

(l0>  Y;  =  gi/KWp    =  A1  e^i-E  , 


A  cosft  coSoJf    B  sin<x  -B  cosa  cos«<+  A  sin  c< 

where  B'  =  -■■  ■ " -— ~ t   and  A'  -  nn    -= f— y— — ;T~"T  * 

K2sinfi  (1-sin^/j  cos2*  )*  '   K  sinrj  (1-ein  |3  com  )  = 

Rt  should  be  noted  th*t  for  incidence  perpendicular  to  the  axis  of  the  boss, 
8  =7T/2  so  that  B*  -  B/K2  end  A'  =/»  A/K2.  For  this  cece,  the  primed  quantities 
(the  TV  mode)  correspond  to  a  wave  polarized  parallel  to  the  s-axis  while  the 
doub]y  primed  quantities  (the  TE  mode)  correspond  to  a  wave  polerized  perpendi- 
cular to  the  z-axis.  Even  for  3  \  ^"/2,  the  problem  could  of  course  be 
formulated  initially  in  terms  of  polarization  components  parallel  and  perpendicular 
to  z  (as  is  usually  done  for  problems  involving  cylinders),  but  this  involves 
some  loss  of  generality  (for  the  reflection  problem  of  the  boss  on  the  plane) 
as  compared  to  the  adopted  procedure  of  referring  the  polarization  to  the 
plane  of  incidence .} 

We  are  now  in  a  position  to  solve  the  problem  of  Fig.  lb  subject  to 
the  second  condition  of  (U)j   i.e.,  we  consider  the  original  wave  and  its 
imsge  simultaneously  incident  on  the  cylinder  end  seek  to  determine  the 
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scattered  field  Ec,    Hc,   by  first  deriving  the  corresponding  if  .     To  facili- 
tete  imnosing  the  boundary  condition  at  r  -  8,   we  first  expand  the  potentials 
in  terms  of  cylindrical  functions  around   the  p.xis  of  the  cylinder.      Now 


iKi'E  =  e-iKr  sin  3    cos   (6  -<*  )  +  iKz  cosp_  e-iK»r  cos   (9  -<*)+  ins 

,.  op 

=  eihz     YZ  inrn  Jn(KT)   cos  n   (T-  9+*  )         , 


iKp-E  =  eiKr  sin$cop   (e+*    >+    iKz  COBf    =  eihzOn£nJn(K'r)cos  n(6+o<   ) 

n»o 


n=o 

where  the  Jn  are  Bessel  functions  ^r.d  £=lifn  =  0  and  £  n  -  2  if  n>  0 
Similarly 

e 
Hence 

V  [+1f"   =  B'eihz£     inf  nJn(KT)  $cos  n(7T-  6+^  )   -  cos  n  (6+^   )J  , 
(11) 

V  "+V"  -  A«eihz  7*     in£    J   (K'r){  cos  n(T-  9+1  )  +  cos  n  (e+o<   )? 
i         p  n  n  t  J 

When  these  plane  waves  are  incident  on  the  cylinder     the  corresponding  scattered 
waves  may  be  expressed  in  the  same  form  as   (11)   but  with  Jn  replaced  by  AH  J 
where  the  H_,    the  Hankel  functions  of  the  first  kind,  are  associated  with  out- 
going cylindrical  wr>ves   (because  of  the  suppressed  time   factor  e~   w    ) .     VJe 
therefore  write 

If'   =  ^"   -fV"  +  tf^  =  P»eihZH    incn[Jn+AnHn~Hcos  r(7r_  °  +0<  )   "  C0E  n(e+*)} 
(12) 

Y"  =V"  +  'V"+Y"  =  A,elhZn    in£    [Jn*AnHnlf  °°6  n(T"  9 +*  *  +  COS  n(e+ol$ 


*.      See  pp.   6,7  in  reference  1.      The  problem  of  a  plane  electromrgnetic   rave 

incident  on  a   cylinder  will  be  found   in   Frank-Mises,   Differentialgleichungen 
der  °hysik,   pp.   869-871,   Rosenberg,   N.Y.,   1943. 
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The  Constanta  Aq,  A   are  now  chosen  to  satisfy  the  boundary  condition 
(Bz  =  J^  «=  °^r=a»  fro*   (7)  and  (g)  it  can  be  seen  that  E^  oc'V'' ,  Bg  oc  -y^i  , 
while  Bz  =  0,  BQ  oc  -^-p   ,  80  that  the  restrictions  on  theV's  are  that 
(  T  '  "  ITr"  =  °^r=a*  ImP08in£  these  °a  (12)  (i.e.,  replacing  r  by  a  in  V"  '  and 
in  ^    and  then  equating  the  coefficients  of  the  trigonometric  terms  of  each  to 
zero)  yields 

(13)  A*  -  -  J  /H       ,  A*  =  -  j'/H* 

n      n'  n  n     n'  n       ' 

where  Ka  sinS  =  K'a  is  the  argument  of  the  Bessel  functions. 

Consequently  (12)  and  the  field  derivable  by  means  of  (7)  and  (g)  repre- 
sent the  solution  to  the  image  problem  of  Fig.  lb,  and  it  can  be  seen  that  this 
is  also  the  solution  to  the  boss  problem  in  the  region  )0l  ^  7T/2  satisfying  the 
required  boundary  conditions  at  r  =  a  and  6  =  7T  /2.  The  latter  is  of  course 
obvious  from  (12)  where  it  can  be  seen  that  V  has  the  same  angular  dependence 
as  Vj  +  "V  and  must  therefore  satisfy  the  same  boundary  condition  at  6  =  77"  /2 
since  this  condition  depends  only  on  9   or  derivatives  with  respect  to  8.  To  see 

analytically  that  the  condition  (E  =  E  =  0)     ,  is  satisfied,  we  note  from 

z    r     0  =  /  2 

(7)  and  (8)  that  E^  oc^  •,  Br  oc  -^~-  while  ^   =  0,  Er  oc  ~~-  .  Hence  E  and 
Br  oc  sin  n  ("7T  /2  ♦<*)  sin  n  (77  /2  -  0)  while  E  oc  ^-g  cos  n(7T  /2  ♦o£  )cos  n(  7772-6), 
so  that  all  three  quantities  are  zero  at  ©  =77/2  as  is  required. 

Hence  when  a  plane  electromagnetic  wave  characterized  by  E  as  in  (5)  is 
incident  on  the  boss  on  the  infinite  plane  the  scattered  reflected  Hertz  po- 
tentials are 

^r  "^p  *  ^c  C^p  "  B'elhE  ^IZ^Vn  6in  a(7T/2+<*)8in  n(^/2  -  6) 
(1U) 

^  r  "^p  *^c  =  ^p  *  A'elhZ  2&ia^nAX  C0B  a(^/2*^)cos  n(7T/2  -  e)   # 
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The  total  reflected  electromagnetic  field  it  then 
(15)  jf  =  IP+IC, 


where  by  means  of  eauations  (7)  and  (8) 

(15a)  Ec=E''+rC=(iK  cosf  &£  +  j  pa*  i£>*  +  (1SL22^  -J3&S  -  i^^gf^^^  . 


Making  use  of  (Ik)  we  obtain 


Ec=  — S — Z  C°8^ Sfco«(A  cosflcosot  +  B  sindK  5Z  **"  AnHn  Bln  n(7T/2+«t  )sin  n(tf/2-  6) 

(l-sin2Qcosi()*(        ^  " 

-(  f-B  cos*  coeoC+  A  sin*}  /Kr  einfl  H  £Z  i**"1  A^H^  n  cos  n<jr/2+oL  )sin  n(T/2-  Ojr* 
L5b) 

+  f  coto([A  cosftcos«<+  B  sin*}  /Kt)aYZ  **'     ARHn  n  sin  n  (T/2  +^  )   cos  n(7T/2-  6) 

+  (-B  cospcosJH  A  sin*  )2  XZ  ^"'''SnVn  C0B  n(T/2  +o(  )   c°8  "("fl"/2-  9M  ^ 
-[sina(A  cos^cow+B  sin*  K  ZZ  iB  AnHn  sin  n{TT/2+<*  )    sin  n(T/2-  ejj^V 


Similarly. 


where 


mm  —  — .  ^ 


(15c)      Hc  =H«C+ 


rc  s(_  ^m+K  cosp^+ (1^^2|%K28ln^ . 


-? 


-19- 

MaMng  use  of  (lU)  again  ve  obtain 
iKz  cosp 

11         7~^ T~l 

'    (1-Bin^  cosV  Y 

X   \  I  (CA  cosflCosw+  B  sin*}  Ar  siiuK  ZZ  *       ARH  n  sin  n(j/2+o<)cos  n(7T/2-  9) 

(15d)  n-l         «   t  "I 

-cosfl(-B  cosp  cos<*+ A  sin*  )2  XL  *  ~  £     A  H     cos  n(ir/2+^)cos  n(JT/2-  6)  I 

+  [(A  cos^coso<+  B  sinrf  K  XZ  i       A'Hn  sin  n(T/2+«t) sin  n(7T/2-  0) 

-  cot^([-B  coe^cosd+A  sinot]/Kr)A  XT  ^  A"H     n  cos  n(7T/2  +  <^)   sin  n(TT/2-  en^ 

■+[siru?(-B  cosficoso(+A  sinol)2  *£Z-  ^i    A  H     co8  n(^/  2  tot) cos  n(7T/2-  9)|ic  > 

where  Ep  and  Hp  are  given  by  equation   (6)   and  where  j[  "  goes  from  1  to  co   while 
2__£      goes  from  0  to  oo.     The  above  constitutes  the  solution  to  the  problem 
for  the  completely  general  case. 

For  K'r»l  and  y}    n,  we  can  substitute  the  asymptotic   form  of  the 
Rankel  functions  into  the  above  series,   i.e.,        H£^(2/i7TK'r)2   i~n  eiK'r, 
H^(2/i7rK'r)2   i"n  +  1  eiK?r.if  K»a<l,   the  series  are  rapidly  convergent 
and  may  be  approximated  by  their  first  few  terms  because  of  the  nature  of  the 
dependence  of  the  coefficients  A^,   a"   on  K'a,   i.e.,   A£  <s  iTT/2  ln(yK'a/2), 
A[*  -iT(K'a)2A,   A^«    -iTn(Ka./?)2tl/(n!)2;      A^s    -i7T(K'a)2A, 

A^-*  i7T(K«8)   A,     A^c-  i7Tn(K«a/ 2)^/(11: J^,      (See  reference  1,  Appendix  II). 

In  a'.V  ■  1.781  . 

2      1/2 
Hence  to  the  order  of  (K'a)  /(K'r)     the  scattered  quantities  may  be 

approximated  by 
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°f'c  ~-  B-    e1Kz   COE<?     (7Ti/2K'r)^    (K'a)2   eiK'r  2  cosmos  0  , 

(16) 

o^"   =  .  A,    e1Kz   cos^    (7Tl/2K'r)*   (K'a)2   eiK'r   (1+2   sin*sin  9) 
'  c  ' 

and 

i(Kz  cos(?+  K'r)  i  , 

°EC  =  -s _ _ —     (tj-  i/2K'r)3(K»a)^  cosfitAcosficoBatf  Bsino02co6cUcoBe  ? 

(1-sin"  p  cos2*)?  C         r  V 

+  (-Bcosfl  coso(-t  Asino0(l+2sin*sin  9)   9 
(17)  P 

-  sing(Aco6flco8<*-f  B6ine<)2co8c<cos6     k    i 

°HC  =  -ylcos^-E^^o-Bl/Bir^)^-  sin^0F*    £j 

where  the  left  superscript  0  hps  been  introduced  to  indicate  that  this  is  the 
approximate  form  for  a  boss  located  at  the  origin  (or  rather  along  the  z-axis) 
of  the  reference  coordinates  (to  facilitate  the  later  extension  of  these  re- 
sults to  distributions  of  bosses  where  the  0  will  be  replaced  by  a  dis- 
crete index  (c)  or  a  continuous  variable  (y)  and  the  expressions  multiplied  by 
a  suitable  phase  factor  chosen  so  that  the  phase  at  the  reference  origin  of  the 
reference  coordinates  is  zero). 

It  can  >e  seen  that  the  formulation  of  the  problem  and  the  derived  ex- 
pressions (5, 6, 10, Li, 15, 16  and  17)  are  considerably  simplified  vhen  the  wave  is 
incident  perpendicular  to  the  axis  of  the  boss,  i.e.,  whenfi=7T/2  and  the  xy-plane  le  the 
plane  of  incidence.  Then  the  component  AF*  (parallel  to  the  plane  of  incidence) 
is  also  parallel  to  the  axis  of  the  boss  and  the  expressions  are  "separable  In 
the  polarizations" — i.e.,  the  previous  priced  or  TM  cuantities  depend  only  on  the 
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B  polarization  component  of  the  incident  wave  which  is  now  parallel  to  k, 
the  vectorial  direction  of  the  Hertz  potential,  while  the  double  primed  or  TE 
quantities  depend  only  on  the  A  polarization  componert  which  is  now  perpen- 
dicular  to  k. 

For  this  case  we  obtain  (on  substituting  fl= 7T/2  in  the  previous  ex- 
pressions), the  plane  waves 

(5.)     e1  =  {-A[sin  (6  -<*)?+ cos(e  -<*)5]+  Bk]e"iKr  cos  (9  -*>  1 

H1  =^JB[sin  (6  -^)r+co3(9^)e]+^e-iKr  cos  (9"*> 

(6')  EP  =     {Afsin  (9  +  c<  )?+cos(6^   )ej    -Bkje   lKr  C0S   (9+°°       ^ 

HP=7^B[sin  (9+o<  )r+cos(9+o<  )e)  i-  Ak"^  eiKr  cos    (9  +  *   >    , 

(10-)      V;  =+(B/K2)   e"iKr  COS   <9  ~*>      ,   V;  =  -   (B/K2)   e1Kr  cos   <9  +  ^  >       , 
If  I  =^(A/K2)   e"iKr  COS   (a  -°<)      ,  V- J  =  /J  (A/K2)   eiKr  cos   (9  +  *  >      ; 

the  reflected  potentials 

(U.)  t>tp^;  =  -  (B/K2)|  eiKr  cofl(9^^  C  ^X  sin  n(7T/2^  )sin  n(7T/2-  9)^ 

Tytytl  =*/U/*i){eiK™B^*k  2  II  iV^coB  n(rr/2+*)cos  n(7T/2-  6)^ 

(which- except  for  the  multiplying  factors, are  identical  with  the  velocity 

potentials  derived  for  the  analogous  icoustical  problems  for  the  cases  of 

0  and  00  surface  impedance — see  reference  1,  page  3)*  and  the  reflected 
field 
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(15')  f  =   Ep+EC        , 


Ec  =   (E^^  +  Eq°  9)  -H^k 

=  -(A/Kr)   A  T~~  I""     A"   H     n  cos  n(7T /2-hrf)    sin  n  (7T/2-  9)   r* 
+  A  2  XZ  in-1£  nAn  Ki  cos  n(T/2  +  °0   cos  n(7T/2-  6)  'e 
-  B  4.  £H  in  A^  H„  sin  n(7T/2+o(.  )   sin  n(7T/2-  9)    'k'  ; 


while 


Hr  =  HPf  Hc  , 

Hc  -   (K'C'>+H,C/9N)+  H*Ck 

=    ^(B/Kr)  A  ZZ  i11"1  An  Hn  n  sln  n(T/2-hoL  )   cos  n  (7T/2-  0)  'r 
-f-^B  ^  £Z  in_1  A^  H^   sin  n(7T/2+o<.  )    sin  n   (7T/2-  9)  'fr 

+mk  2  £"  in^n  An  ^n  C0F  "C^/24  *  )    cos  n   (TT/2"  e)    ^ 
For  Kr»  1,  Ka<,l,   the   scattered  Quantities  reduce  to 
(16')  °VC   =   (B/K2)    (7Ti/2Kr)^   (Ka)2   eiKr  7  cosfX  cos  9  , 

0Y^  =  -/^  (A/K2)(7Ti/2Kr)^   (Ka)2  eiKr  (it  2  sin*  sin  9)    j 
(17')  °EC  =  -  eiKr(7T  i/2Kr)^(Ka)2[A(l+2  sino<.sin  0)%  -  B  2  cos*  cos  9^] 

(Expressions   (17')   are  essentially  those  obtained  by  Rayleigh6   for  the  above 
problem.) 
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On  adding  the  p  and  c  terms  we  obtain 
(17")  f   =  e1Kr  C°G  (9M  }{a  sin  (9 +*  )  ?ta[cos(9+o<  )-qQ  elq]  ^  -  B[l-Ok  ela]  ^ 

Hr  =  *>  eiKr  C0S  (6+^)^B  sin(9+^^[B  cos(9+oO-Qk  eiqJ  3  f  AJ>0Q  eiqJk^ 

where     q  =  Kr  [l-cos(9+<<  )]  +  77"  /u  -   2Kr  sin2  [(9+^  )/2]  -r  7T/4, 
0  =  (7T/2Kr)2  (Ka)2|  (1+2  sin*  sin  9)+  tf'(Ka)2  ...  J    , 
Q^  -   (7T/2Kr)?  (Ka)2}  2  cosotcos  9+  &  (Ka)   ...    ?  . 

(The  notation  (•(Ka),  ^(Ka)  ,  merely  indicate  that  the  next  order  terms  in 

2 
the  curly  brackets  are  to  the  order  of  (Ka)  or  th?.t  the  next  order  terms  in 

Q  are  to  the  order  of  (Ka)'*). 

Further  consideration  of  the  implications  of  the  above  will  be  reserved 

for  the  last  section  of  this  report. 

B.  Distributions  of  Bosses. 

The  asymptotic  solutions  obtained  for  the  single  boss  can  be  extended 
to  yield  approximate  solutions  to  some  problems  of  scattered  reflection  from 
certain  planar  distributions  of  such  bosses.  The  initial  -ssumption  is  that 
the  spacing  of  any  two  bosses  is  greater  than  the  incident  wavelength  so  th?t 
the  excitation  of  each  boss  is  essentially  only  the  incident  plane  wave:  i.e., 
the  contributions  to  the  excitation  of  each  boss  arising  from  the  waves  scat- 
tered by  the  remaining  bosses  is  neglected — the  so  called  "single  scattering" 
hypothesis.   (This  will  be  further  discussed  in  section  5-) 
l)  Patterned  Distribution. 

Consider  a  wave  characterized  by  E  (as  in  (5))  incident  on  2N-V-1  =99 
bosses  of  radius-a^ whose  axes  lie  in  the  yz-plane  parallel  to  the  z-axis  and 
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are  spaced  symmetrically  along  the  y-axis  with  spacing-b— a  reflection  grating 

of  bosses.  Since  the  procedure  for  solving  this  problem  is  essentially  the 

same  as  that  of  section IIB  reference  1,  it  will  be  but  briefly  sketched  here." 

""e  consider  the  point  of  intersection  of  the  axis  of  each  boss  with 

the  y-axis  of  our  reference  frame  to  be  the  origin  (located  at  0,  qb,  9  in 

the  xyz-reference  frame  of  the  central  boas,  q  =  0)  of  a  cylindrical  coordinate 

system  r  ,  9q,  z,  in  which  the  distance  of  a  field  point  P  may  be  expressed 

as  R  =  r  (cos  0  i  +  sin  9  'J)  -t  zk  =  r  cos  9  i+ r  sin  9  j-+zk.  The  Hertz 
-q    q       q        q  q      q 

potentials  of  the  incident  and  image  waves  for  the  q*th  boss  differ  from  those 

for  the  central  boss  (as  in  (10))  only  in  that  R  is  replaced  by  R  and  the 

-iK'qbsin^ , 
expressions  multiplied  by  a  phase  factor  e  (as  obtained  by  setting 

8  =  7772,  r  =»  y  =  qb,  z  =  0,  in  either  e  -i'-  or  e  -P*-)  to  correspond  to  the 

Dhs.se  being  zero  at  the  central  boss.  Similarly  the  scattered  potentials  (the 

■^c   of (14.)  or  (16))  differ  only  in  possessing  r  and  9  (as  will  be  indi- 

q 

cated  by  writing  if  )   and  in  being  multiplied  by  the  phase  factor.   Hence 

the  scattered  contributions  of  the  entire  distribution  may  be  expressed  as 

q„>   -iK'qb  sin* 
e 


as)  y.'iz  „  V. 

a  =  -  N      c 

Now  if  the  extent  of  the  distribution  is  «  r,(i.e./?7  b«r)then  we  have 

approximately  r  «.  r-qb  sin«<  ,  9q-%9  -  (qb/r  )  cos  9  -}  9  .  On  embodying 

q  • 

these  approximations  into  the   y"  and  substituting  the  asymptotic  repre- 
sentation for  the  Hankel  functions  (neglecting  qb  sin  9<<  r  in  the  denomin- 
ator) we  find  that  the  resulting  sum  over  q  .is  a  geometrical  progression  which 
can  be  summed  explicitly.   If  the  restriction  K'a<l  is  then  invoked  we 
obtain 

9.  The  procedure  was  originally  applied  to  the  problem  of  a  plane  wave 

normally  incident  on  a  grating  of  parallel  cylinders  by  C.  Schaefer  and 
?.  Reiche,  Ann.  d.  Phys.,  21  817  (1911). 
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(19)  ^c  =  °Yo  sinW/Vsin  P        ,  T    =  K»b(   sin  6+ sin*  )/2      , 

where  the  °Vl  are  as  In  (16)   or  (16').     We  therefore  obtain  for  the  scattered 

'  c 

reflected  field  components 

(20)  f  =  ZP+°f     sinW/sinT      , 

where  F  »  E  or  H  and  the  appropriate  quantities  may  be  found  in  (6)  and  (17) 
or  (6')  and  (17»).  (That  (19)  implies  (20)  can  be  immediately  verified  from 
(15a, c)  by  noting  that  for  K'r  »  1  the  F°   depend  essentially  only  on  the 
1^%e   and  their  derivatives  with  respect  to  r -•  the  derivatives  with  respect 
to  9  which  vanish  as  v~^'      being  neglected.  The  quantities  (20)  could  of 
course  also  be  obtained  by  operating  directly  on  the  °FC  as  in  (18)-(19).) 

The  above  may  be  readily  extended  to  a  superposition  of  several  such 
reflection  gratings — which  are  parallel  and  coplaner  and  overlap  only  at  the 
central  boss  and  differ  in  the  number  of  elements  {$[) ) ,  their  spacing  (b) , 
and  their  radii  (a) — by  summing  the  contributions  of  the  various  gratings. 
We  would  then  obtain  (20)  but  containing 

(2D        Hj  *j  f  (sin^j^/.inPj)  -y^l 

where  'y  ^   -  0  for  one  particular  grating   (that  of  largest  a*)   and  47 ,  -  1 

for  all  the  others  to  condensate  for  ih*  overlap  at  tho  cantor,     M  tho  specular 

ancle  of  reflection  (essentially  whoa  C  «  -  <X)  tho  W     «  0  for  all  j  so  that 

si*'??,  nJni* r  j  -*%  ••*  (a)  3ri»io«  EZj  *j  ^ r  y \ h whleh  for  m11 

a,  «  a  yields  a  £Z  (9?  j-  y  J  ■  *  ^7, where  *?7  is  tho  total  ncusber  of  bosses 
la  the  superposed  distribution.  We  therefor*  obtain 

(22)  (yr  =  3**  V«L  . 

—    —     —     tr*  -  on  , 

This  procedure  nay  also  be  extended  to  consider  overlaps  occurring 
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else*here  than  at  the  center,  i.e.,  when  the  spacing?  of  various  gratings  are 
cor.mensurate,  by  compensating  for  the  overlaps,  es  above,  by  the  addition  and 
the  subtraction  of  the  missing  terms.  It  can  be  seen  that  "t  the  specular 
angle  these  terns  are  equal  to  1  so  that  re  again  obtain  (22).  Also  since 
this  is  true  for  any  distribution  which  may  be  decomposed  into  patterned  distri- 
butions, we  may  therefore  conclude  that  this  is  also  true  for  any  distribution 
of  parallel  bosses  of  equal  radii— subject  to  the  initial  restrictions  and  to 
the  same  order  of  approximation  (viz.,  the  single  scattering  hypothesis  and  that 
the  extent  of  the  distribution  be«  r)—  i.e.,  that  at  the  specular  angle  the 
effect  of  such  a  distribution  may  be  found  by  multiplying  the  effect  of  the 
central  boss  (°Y'  )  by  the  total  number  of  bosses. 

2)  Uniform  Random  Distribution, 

Consider  F  incident  on  a  uniform  random  distribution  of  parallel 
bosses  of  equal  radii  on  an  infinite  yz-plane;  the  axes  of  the  bosses  being 
distributed  along  the  y-axi s  according  to  a  density  function  M  which  is  inde- 
pendent of  y,  so  that  M  dy  equals  the  number  of  bosses  in  area  of  width  dy  by 
unit  length  (or  ersentislly  in  dy) .   Analogous  to  the  summation  procedure 
of  the  preceding  section  and  based  p.e   previously  on  the  single  scattering  hy- 
pothesis we  must  nor  evaluate 

(23)         -v„=  ^r^"'7^"* 

where  the  superscript  y  indicates  that  in  (16)  for  ° f „   we  replace  r,  e  by 
r',  9'  which  pre  functions  of  y — i.e.,  the  running  coordinates  from  a  point 
on  the  distribution  to  a  fixed  field  point  P  (while  r,  G,  will  be  retained  to 
specify  the  position  of  P  in  the  reference  coordinate  frame  at  the  center 
of  the  distribution). 
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(a)  Small  Finite  Distribution  . 

For  a  finite  distribution  extending  from  y  =  -  d/2  to  y  =  +  d/2  (d«  r) 
we  obtain  by  essentially  the  same  procedure  as  for  the  grating  (i.e.,  replacing 
r'  by  r  -  y  sin  6  and  9'  by  6,  neglecting  y  sin  9«  r  in  the  denominator,  and 
evaluating  the  resulting  simple  integral  as  in  reference  1,11  CI)  that 

(2u)  f;  =  °%  Md  sin  p/r 

where  f  is  as  in  (19)  with  b  replaced  by  d.  Hence 

(25)  Fr  =  fPf0FC  Md  sin  T/r 

This  may  also  be  extended,  by  summation,  to  a  superposition  of  such 
distributions  which  are  parallel  and  coplanar  and  where  each  is  characterized 
by  o  different  density  function,  length  and  boss  radius.  Then  provided  that 
all  aj  =  a  we  would  again  obtain  expression  (22)  at  the  speculer  angle 
(fj  =  0,   sinPj/f1  j-»  1)  wlth<>?=  Md  or*??  =  £Z  Mjdj  for  the  total  number  of 
bosses  in  the  distribution. 

(b)  Infinite  Distribution  . 

For  a  distribution  extending  from  y  =  -  oo  to  y  =*oo  we  evaluate  (27) 
at  an  arbitrary  constant  field  point  P,  whose  coordinates  in  the  reference  frame 
are  x  =  D  =  r  cos  9,  y  =  C  =  r  sin  Q,  z  *  z,  by  the  method  of  stationary  phase 
(section  IIC2,  reference  1).  From  (23)  we  obtain  on  substituting  from  (16) 

(23')  V'c   =  B'elhz(7T  I/O.)}  (K.a)2  2  cos^f00  .«,^,-7  ^  >  (r.)-^coseQ  dy 
y.  =  .  A..ih*(Tri/2K.)*(K«e)2Iir00  J*'^'*   si»*  >(r')4[l+2  sin.  sinewy 
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Th  e  integrals  may  be  expressed  in  terms  of  8*  as 

A/2 


(26)      D*  e-1K,c  *»+y/2   eiK'D(sec  e«+  tan  e*  sin*)  £    -J  jA  e,  d  Q, 


(where  the  bracket  is  either  that  of  (23»)  or  (23"))  and  yield  on  evaluation 
(the  phase  being  stationary  «t  9'  =  -oO 

<2?>  r  ]e.  =  -* eiK,r  cos(9+* )  -^(2iff A*)* .  w^t/8)   . 

Hence 

(28)  'V''   =  IB'   e1£P*-  K'M77a2  2  cos2<X     sec  *     , 

c 

V*"  =  -1A»   e1"?"-  K'M7ra2   co»2*  sec«(       1 
c 

and  the  scattered  reflected  potentials  are 

(29) 

f*  -Tf"^"  =-V-"    (1-1K-V2  cos  2oCseco<)    , 


where  V  =  M7Ta2/2  andV' ',  Y'"'   are  RE  in   (10)*     The  fleld  as  °btalned  f«« 


P'    '   P 
(15a, c)   is  then 
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¥  =  IPH-  EC     , 


Ec  =  - 


__J2K!lJLZ!ii_^-  i  [  cos«  (Acosq  cosd  +  Bsine*  )cob(0  +* )2cos2*< 
cos<X  (l-sin2|3  cos2*  )H  ^  * 

■+  (-Bcob/3  cos^-f  Asin  c<)sin(e+ol  )cos2<T]^ 

2 
-t[-cosfl  (Acos*  co8^Bsin<X)sin(e-h/)2  cos<* 

^  (-Bcoag  coscH  Asin*)  cos  (9+o<  )c©62o<Je 

-sing  (AcoeG  costf+  Bsirio(  )2cos  o(  k  J 
(30)  I 

=  _  i2K«V  seco<[    2  cos2*  E«P+cos2o<E"P] 


* 

9 


Hc  _  .  .^i2^^^--— -r([Uc<>s3  COBJ+  Bsin*  )sin(9^  )2cos2* 
"  cosoCO^infc    cos2*  )HL  1 

-cosrt  (-Bcosa  cos**  Asin  <*)  cos  (6+o<  )cos2«[]  r 
-H(AcosA?cosfl<+Bsin<*)co8(e-*<02co6  * 

+   cos/3  (-Bcosp  cosJ+  Asin^)sin(e+ot)coB2^J 
+  sin/3  (-Bcosq  cobc(+  Asin  <ji)cos2^.  k  > 

=  -i2K'V  sec*[  2cos2<*  H,p4  cos2<*H"P3       j 
where  Ep,   Hp,   are  as  in   (6). 
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For  the  case@  ^TT/2   (incidence  perpendicular  to  the  axis  of  the  boss) 
the  above  simplify  considerably.     Thus 

r  -  -    (B/K  )   e  '    (1-1KV4  cos*)  , 

(29') 

V  ;  --*?  (AA2)    ei*r  cos   (e«  )    ^^  co8  2Ueeco() 

(which  except  for  the  multiplying  constants  are  identical  with  the  velocity 
potentials  for  the  cases  of  0  and  oo  impedance— p.  17,  reference  D;  and 
the  field 

jf  =  EP+  Ec 

c-c  _       mOVV  ^iKr  cos   (8M)         ,(  r  ^  >o  ~>  -> 

£    "  "  l2KV  e  ;£eco<^Acos2^sin(e+^^^co8(e^)g-B?cos2^   ^k? 

=  -12KV  sec<*rCos2o<  eJ+  2  cos2*   E?   1 
(30')  L  "n  _J-  J  * 

f  =  jf 4  iiC  , 

Hc  =  -  i?2KV  eiKr  C0S^V  )6ec^B2co6^   ^(^ft,  cob(0**  >3>eo82^ 
=  -  12KV  secocl"  2  coe2^  HP+  cos2ol  Hj  1 

where  EP  and  HP  are  as   in   (6»)   and  "here  the   subscripts  |J,jL,    indicate   the 
components  fl  and X to  the  plane  cf  incidence.    (Note  the  relationship  of  the  c 
and  p  nuantities  of  these  expressions,   from  which   it  can  be   seen   that  were 
either  A  or  B  eoual  to  zero  then  the  scattered  component   (c)   woulr   be  directly 
proportional  to  the  specular  component   (p) . 
On  adding  the  p  and  c  terms  re  obtain 

E1"  =  [A(l-i2KVsec<Cos2o<  )[sin(e+rf  )rfcos(e^  )eJ-B(l-HKV  cos*  )^iKr  cos(e  +  *) 

Hr  =  ^B(l-HKVcos<<  )[ sin(9+*  )r+cos(e+<X  )e]+A(l-i2KVsec*  cos2*  fk?eiKr  coe(9*<0 
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There  is  little  resemblance  between  the  scattered  components  of  these 
expressions  end  those  of  the  finite  distributions  (JO),  (25)  treated  pre- 
viously. The  chamct.ri.tic  ten*  of  the  single  semicylindrical  boss  no 
longer  appears;  in  feet,  -  have  lost  .11  trace  of  cylindrical  vaves  and  the 
.nattered  reflected  .eve  consists  of  .  superposition  of  plane  .eves  of  dif- 
ferent amplitudes  end  phase  traveling  in  the  specular  direction.  The  symbol 
V,  the  volumetric  departure  from  the  plane,  is  introduced  to  eliminate  ell 
trace  of  the  boss  from  the  expressions.  Since  TT.2/2  is  the  volume  of  a 
semicyllnder  of  unit  length,  V  represents  the  total  volume  of  bosses  per 
square  centimeter  of  distribution  and  has  the  dimension  of  length. 

For  a  superposition  cf  such  distributions,  ehich  are  all  parallel 
and  coplanar,  me  obtain  the  above  expressions  containing  By     «*"«>  de" 

notlngE  V  «-  ™*  "1»'«10  d'P"rttt"  fr°°  °"  Pl0n''  ^  '  yleld'  ""  '"" 
expression,  as  for  the  single  distribution. 

».  could  of  course  have  obtained  (30)  by  performing  the  above  operations 
(26),  (28)  directly  on  the  field  quantities  as  given  by  (17).  The  resulting 
expressions  mould  then  be  In  term,  of  the  unit  vectors  associated  mith  the 
stationary  value  rather  than  in  terms  of  the  reference  unit  vectors  and  .ould 
require  transformation  «,   facilitate  comparison  vitn  the  plane  save  compon- 
ents. We  mill  demonstrate  the  procedure  for  future  reference. 
We  wish  to  evaluate 

C°°     v  c  -iK'y  sin*  , 
(3D  EC  =  «J     ^  e  dy    . 

For  greater  lucidity  we  will  first  perform  the  trans for™ tion  and  then  the 
integration.  For  brevity  we  write 

(32)  ?EC  =  Ejt'+^Sh  E$c      , 
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the  scalar  coefficients  «,,  %,   Z}   („hose  ralues  ca„  „.  ^  ^  ^^ 

-Hh  (17))   being  functions  of  the  variables  r-   and  ...     „  m  transforn, 

this  entity  ana  express  it  in  te™  of  the  te,e  unit  vector,  of  .be  reference 


frame;   r,   9,  k.     Thus 


7*c  -  y„c 


*s  may  be  seen  from  the  sketch  below,  the  scalar  products  are  equal  to 

m         tf*  -  cos  (9  -  9.)    ,    S..9  a  C0B{yT/2   _  (Q  _  effl  m   sin  (q  _  ef)  ^ 

r'.^=coS[7T/2-K9-e.)]=-  sln(6-9.)   ,  9-.$  -  cos  (9  -  *•) 


? 


4*  5 


>r 
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Consequently 

(35)  yEc  =  [f^cob^  -  e«)4E2Bin(e  -  9'  jr+^sMe  -  e»)+R2cos(8  -  e«)I  9+E^ 


Expression  (31)  is  then  equal  to 

ihz 


Ec  = S 


j-  frl/2K')i(K'*y2*j"    eiK'(r'-7  sincO(rf)-l/2 


(l-sin2^  cos2*^  )*  -/-oo 

(Acosq   cosof+  B  sino(  )2  coso<  cos  9  cos    (6  -  9") 


J   [c°8^ 


V> 


1 

+  (-B  cosScoerf+A  sin*)  (1+2  sin*  sin  6»)  sin  (9  -  9»)1  "r 

+  J-cosR  (A  cosft  cos*  +  B  sin<*  )2  cos*  cob  9f  sin  (9  -  9*) 
+  (-B  cos/*  cob^+A  6in*)(l+  2  sin  *  sin  9')  cob  (9  -  9»)J  3 
-  sing  (A  cosft  COBW+  B  sin*  )2  cos*  cos  9'  k£  dy 
which  on  evaluation's  previously,  yields  (30). 

C.  Distributions  of  Cylinders  . 

The  results  of  the  preceding  section  can  be  readily  extended  to  the 
problems  of  E^  (as  in  (5))  incident  on  similar  distributions  of  whole  cylinders. 
The  expressions  for  the^  and  FC  for  the  small  finite  distributions  (as  in 
(l9)-(25))  are  formally  identical  with  those  that  would  be  obtained  for  the 
analogous  distributions  of  cylinders  if  it  is  understood  that  °j'    and  °FC 
represent  the  approximate  solutions  for  the  single  cylinder  at  the  center  of 
the  distribution  for  K'r>>  1,  K'*<1. 

The  solutions  for  the  single  cylinder  may  be  readily  obtained  by  what 
is  essentially  a  simplification  of  the  previous  procedure.  Thus  from  (10) 
we  obtain  the  forms  of  the  j.   and  consider  them  to  be  incident  on  the  cylinder. 
The  total  potentials  are  then 
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T  --t±  +  t  c  =  C  eihz  Yl  SB  Q[jn+  ^  *Q   cosn(T-e  +  <)     J 

where  f or  W ,  C  =  B',  A  =  A'J  while  for  V".  C  =  A1,  A  =  a".   The  scattered 
'  '      *  n    n  •  *      '  n    n 

potentials,  7"  ,  are  then  eoual  to 

V"  ;  ■  B.  elhZ  H  iD£n  A;  Hn  cos  n  (7T-  0«)  % 

Y"'  =  A'  eihz  5~"  ine  „  A"  H  cos  n  (ff-  6  +  JL ) 

c         ' —    n  n  n  ' 

•   " 
where  the  A  ,  A  ,  are  as  in  (13).  The  rest  follows  as  previously  and  to  the 
n   n 

order  of  (K'a)2/(K'r)2  we  obtain 

°f   =  B»  eihZ  +  iK'r  (i7T/2K'r)^  Fl/ln(  YK'a/2)  +  (K«a)2  cos  (9  -<<)] 
'   c  > 

(37) 
of"  =  _a«  eih2  +  1K'r  (i7T/2K'r)^  (K«a)2[l+2  cos  (6  -  «)1 /2     J 

and  the  field  as  obtained  from  (15a, c) 

ihz  +  iK'r 
e ^__ 

(l-sin2fi  cos2<*  )2 

x.Scose(Acos<?cos«y+Bsino()[l/ln(rK'a/2H(K,a)   cos(9  -<*>]? 
+  (-Bcos^  coM  +  Asino<)(K'a)2([l+2co8(e  -at  )}  /2) 'q 
-  sinfl(Acos£cos<rf  +  Bsino<)[l/ln(>'K'a/2)  +  (K'a)2cos(e  -<*f|1c? 
°HC  =  -  ??f  cosfl  °Ece/r^  (°E^/sin^  )  'S  -  sinp  o^g^? 


ihz  +  iK'r  , 

V» S_ _ .  (i  7T/2K'r)5 
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Hfld  the  incident  wave  been  E^  (as  in  (6))  the  above  would  be  modified 
only  to  the  extent  that  cos(9  -o<)  would  be  replaced  by  -co6(9"t-o()  and  that  the 
bracket  containing  the  In  term  (that  of  if*  )  would  be  multiplied  by  -1.   (It 
should  be  noted  that  the  (K'a)  term  in  the  bracket  containing  the  logarithm 
is  included  only  to  facilitate  comparison  with  the  previously  derived  ex- 
pressions for  the  bosses  and  that  its  contribution  is  negligible  compared  to 
the  logarithmic  term).  It  can  then  readily  be  seen  th~t  if  these  two  cases 
are  added  (i.e.,  E  and  T    simultaneously  incident  on  the  cylinder,  the  identical 
image  problem  treated  for  the  bosses)  the  In  terms  cancel  while  we  have 

cos(6  -oO*cos(9-H.)  =  2  cos  6  eosot  for 'V'c  .and  -e  have  for'V  "   that 

c 

£l4  2  cos(6  -oO  +  l  -  2  cos  (6-K)l  /2  =  H-2  sin  9  sin  <K   ,  as  was  of  course 
expected.   (See  expressions  (16)). 

For  the  infinite  distribution  and  an  incident  wave  J~,  we  obtain  by  the 
method  of  stationary  phase  applied  to  (23)  utilizing  the  potentials  of  (37),  that 

T  ^  =  iB«  e  ^  ~  K»  M7Ta2  sec  ^  Q»   ^ 
(39) 

7  „  =  i  A'  e  -q  ~  K»  M7Ta2  sec<*  0   • 

where  for  the  reflected  wave 

(40)      q  =  p,  Q»  =  l/(K'a)2  ln(  y^'a/2)i- cob2o(  ,  Q"  =  [  l+  2  cos  2o(l  /2 

while  for  the  tra.nsmitted  wave  (obtained  by  replacing  9'  by  -  9'  before  evalu- 
ating the  integral) 

(U)         q  =  i,  Q.»  =[l/(K'a)2  ln(jTK'a/2)  -  l  I    ,  Q"  =  -  l/2 
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Had  the  incident  wove  been  Ep  rather  than  IT  the  sbove  would  be 
altered  only  to  the  extent  of  interchanging  i  and  p  in  q  and  multiplying 
0'  by  -1.  As  previously  the  addition  of  these  two  cases  may  be  used  as  a 
check  on  the  results  derived  for  the  bosses.  Thus  for  x)0   (q  =  p,  the  re- 
flected wave  for  E  and  the  transmitted  wave  for  Ep)  the  In  terms  cancel  to 
yield  cos  2<X+1  =  2  cos2^  for  V*c  while [(1  +  2  cos  2o(.)  /2]-  l/2  =  cos  2c( 

forV"  as  obtained  previously.   (See  expressions  (28)). 
c 

The  corresponding  field  quantities  may  be  readily  obtained  from  the  above  by 
differentiation  as  in  (15a, c) ,  or  by  inspection  of  (30).  Thue^  FC=iK'Mfa2(0l?t\Q"F"<5) 

(It  should  again  be  noted  that,  in  the  In  bracket,  it  is  the  In  term  that 
predominates  and  the  other  term  which  is  negligible  ie  included  only  to  facilitate 

comparison.  It  should  also  be  noted  th?t  M7Ta2  is  the  volume  density  of  cylinders 
or  essentially  the  volume  of  material  per  centimeter  of  distribution  along  the 
y-axis) . 

U.     The  Bern 1 s pher i cal  .Boss  . 
A.  The  Single  Boss  . 

Consider  a  plane  electromagnetic  plane  wave  characterized  by  E  incident 
on  a  hemispherical  boss  of  radius-a  on  an  infinite  xy-plane  as  in  Fig.  2a.With- 
out  loss  of  generality  the  plane  of  incidence  is  taken  to  be  the  xz-plane.  Fe 
seek  a  solution  of  the  problem  as  formulated  in  section  2  subject  to  the  boundary 
condition  th^.t  the  tangential  component  of  E  vanishes  on  the  plane  and  boss, 
or  that 

U2)       (EY  =  E  =  0,  or  E  =  E  =  0) 

*   r        y    x    6  =7/2,  z  =  0 

(Ey=  E  =  0) 

v  «{   9     r  =  a 
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Fig.   2(a) 


Pig.   2(b) 


H  (region   of  Z  negative) 


I  (region    of    Z   positive) 


Figure  2 

(a)  A  plane  wave  incident  on  a  hemispherical  boss  on  an  infinite  xy-plane. 
The  plane  of  incidence  is  the  xz-plane  in  which  $i  mekes  an  angle  o(  with 
the  z-axis. 

(b)  Two  plane  waves,  $j  and  its  image  in  the  xy-plane,  np,  simultaneously 
incident  on  a  sphere  at  the  origin. 
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The  image  thechnique  discussed  and  employed  previously  will  be  employed  and  we 
will  concern  ourselves  with  the  problem  as  represented  by  Pig.  2b  and  seek  to 
determine  the  Hertz  potentials  of  the  form  II  =  j r. 

In  general  the  incident,  field  can  be  initially  'lecomposed  into  components 
parallel  and  perpendicular  to  the  xz-plane  of  incidence,  expressed  in  terms  of 
the  reference  cartesian  coordinates  and  then  transformed  to  the  spherical  co- 
ordinates of  the  boss,  as  follows  sequentially: 

E1  =  AE.J  +  BE*  =  ^(-cob*!*  sin*^)*  Bj'Je1-1*- 

=  e  -i  ~   $[A(-co8«<.  sin9  cos^+sinc*.  cos9)  -fB  sin9  sin^Jr* 

t[-A(cob«(  co89  eo8<f+sln*  sind)  +  B  cos9  sin<f]  Q 
U3) 

+  Ta  cosc/sinV+B  cos<^"]^f 

and 

VT  -  <*j  n  y  E1   =/^[  B(cos*  i   -  sin*  k)+  A j  $  e  -i  - 

=  /y>  e1-^'-  \  fB(cosot  sin  9  cos  <^   -  sinrf  cos  9)-r  A   sin  9  sin  */] 

-r  £b(cos<*cos  9  cos  <f-r  sin  <<sin  9)    -V  A  cos  9  sin^l   9 
■+£-B  cos*  sin<^A  cost^](^> 

In  these  expression  the  position  vector  of  a  field  point  P  is  r  «  r  (sin  8  cos^i 
t-  Sin  Q  sln/'J  ♦  cos  §  k);   the  incident  propagation  vector  is 

X,   -  In,  -  -  KsinotT  ♦  coi<Ak)  ; 


*? 


^> 
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and  the  transformations 

i  =  sin  0  cos<{  r^  cos  9  cos<{  9  -  sinifi^  ^ 
j  =  sin  9  siny  r+cos  9  siny  9  +  cost{l( 
k  =  cos  9^  -  sin  9  9     , 

have  been  employed.  It  can  be  seen  that  H*  may  be  obtained  from  E  by  re- 
placing A  by  -B  and  B  by  A  and  multiplying  the  resulting  expression  by  yi   . 

The  image  wave  in  the  xy -plane  (characterized  by  EP  and  HP)  is  the  specu- 
larly reflected  wave  that  would  be  obtained  were  the  boss  absent  and  is  deter- 
mined by  imposing  the  first  of  conditions  (U2)   on  E/--frEP.  Thus  on  replacing 
Kj  by  the  image  propagation  vector  K_  =  KnL  =  K(-sin<*  i  4cos*  k)  we  obtain  Ep 
by  changing  the  signs  of  the  i  and  j  components  of  E  ,  and  Hp  by  changing  the 
sign  of"  the  k  component  of  B  .  Hence  we  obtain 


EP  =  JA(cos,/i4sin.('k)   -  Bj£e  _p*- 


iK„-r 


and 


=  e  ~ P  ~  }[A(cos  «<.sin  9  cos<f4  sin  e*.cos  9)   -  B  sin  9  sincfj'r 

-*•  f  A(cos<*cos  9  cosy -sin  <*.  sin  9)   -  B  cos  9  siny]  9 
(U) 

-[A  cos*  siny+B  cos<<]y  r 
Hp  =  ^^x  EP  =/^B(cos/i  -tsin/k)  +  Aj^  e1-^'- 

-yfo  e  ~P  -^[B(cosot  6in  9  cosy+  sin*,  cos  9)+  A  sin  9  siny]r 
+fB(cos«C  cos  9  cosy  -sin  «<sin  9)  +  A  cos  9  sin  yj® 
■+  £-B  cose*.siny4A  cos  <(]((/ 
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To  proceed  further  and  discover  the  field  E8  and  H  scattered  by  the 
sphere  subject  to  the  boundary  condition  at  r  =  a,  we  now  find  it  convenient 
to  introduce  two  Hertz  potentials  of  the  form  II  =  jr.     Substituting  into 
(3)  yields 

V  xVx(fr)  -Vv,'(^r)-^'Vr  =  0 
=7  x(7f  ff  4^)-V(r«T7f  )-V(^'l)-K2tr  =  0 

=  (VV)V  -I-lV2^  (r-V  )?t   -(V^-V  )r-(r.v)V^4  (VV  -V  Jr-rV^xft^) 
"t  £  x  (^^h-(V^)V  -r  -ifWr  -  K2^  "  0  -* 
=  (V2+K2)^=0 

so  that  the'Vare  the  usual  spherical  eigenf unctions.    We  therefore  obtain 
from  (1)  and  (2)  that  the  field  can  be  expressed  as  the  sum  of  the  two  partial 
fields 

E«-  (K2+VV  •)(^'r)=  K^'r+^C^'rJv'  •'r+r-V  (^r)  =  K2^'r-fv[^  ^"^] 

(45) 
H'  =  -  iK/^^xC^r)  =  iK/^fxV^'r) 

and 

E"  =  i^CwVxC^"r)  =  i/uco('e  i-~  1—  -  v?  i  i~?  (^"r) 

(46) 
JT  =  (K^7  -)(t"r)  ^  -r  (*2  +  ^>)^  J  ^  +  ?  ~^  -£^  (f«r)  . 


10.  A  discussion  of  the  functions  and  expansions  used  in  this  section  will  be 
found  in  Stratton  (reference  5)  Ch.  VII. 
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where 


E1  =»  j/1*  E'P+E.'8 


E"  =  E"*-*.  E"P+  E"8 


H'  =  H'^  H'P+  H,s 


v  =^+^P+^; 


H"  =  H"i+  H"P-t  H"6 


will  now  be  determined.  It  should  be  noted  from  (4.5)  and  (46)  that 

Er  =  E^.  =  (K2+D^)(^-«r)  and  Sr  =  H"  =  (K2-*  D2^"*")  (where  Dr  =  T~2  »  for 

brevity^  The  following  procedure  hinges  on.  these  two  expressions  since  the  radial 

components  of  the  incident  and  image  waves  yield  the  corresponding  jr's 

with  which  we  determine  the  scattered^  » s  and  hence  the  scattered  field. 

It  should  also  be  noted  thnt  the  primed  quantities  (of  the  electric  type)  have 

no  component  of  3  in  the  r  direction  and  are  therefore  the  TM  mode,  while 

the  double  primed  quantities  (of  the  magnetic  type)  have  no  component  of  E 

parallel  to  r  and  are  therefore  the  TE  mode — transverse  referring  of  course  to 

the  vectorial  direction  of  II. 

The  procedure  from  this  point  on  is  slightly  more  complicated  than  for  the 
analogous  cylindrical  problem.    We  determine  the  associated  Hertz  potentials 
from  the  radial  components  of  (43)  and  (44) •  Now 

i^i'I     -iKr(sin  9  coswf  sin^cos  8  coso<.) 


00  n  .        . 

=  EI  in(2n+l)Jn(Kr)  JZLtm  f~M  ^MWow^  )cos  ■  *<-!/ 


n=o 


m=o 


(n-r  m)  I 


s  H  nM  (-D1 

«— —•         nm 


11.     The  treatment  of  the  problem  of  a  plane  electromagnetic  wave  incident  on  a 
sphere  on  which  the  essentials  of  the  following  development   is  tesed  will 
be  found  in  either  Frank-Mises   (reference  8)    pp.   871-375  or  M.    Born,   Optik, 
pp.   274-?85,   Edwards,   Mic-iigan,   1943. 


-&- 


( where  the  1  are  the  spherical  Bessel  functions  in  Morse's  notation  and  the  K 
v         Jn  n 

are  the  associated  Legendre  functions)  while 

iKr-T     iKr(-6in  6  cos</  sinot  ■+  cos  6  cosoc  )    <r —  __   ,  -,\m 
e-P-=e  ^  =2—  IInm  (_1'  ' 

(The  symbols  Hrm   above  and  the  host  that  follow  are  introduced  for  brevity 
and  will  be  discarded  as  soon  ^s  they  have  served  their  purpose).  Conseau«ntly 
we  can  formally  express  the  radial  components  as 

r   KA   »4  T  sin<*  W    '  iKr 

=     i^IZi^n+l)jnIZ£ra{SfV[^C0S4   AP£(coso<  )„    cos.* 
n=l  m=o  v  '  t 


B 
sin 


-  P^  (cos^  )   m  sin  m  (^    (-!)«     =    J^     Gn*   <"«"  , 


US) 


iK   *r 
r  (A  JdT*  sin*       Di<  '       iKr    "  "    *—    Gnm  (_1) 


A  }       x   e 


i1  »  07  (_B  ^—  ♦  -r=—    -I ) 

r       /  v        **.       sin  ok      3</  ' 


iK  -r 

1   - 


Or 


=  ^^in(^n+l)    j^f-fefiyt  P>o8  9)j    -BP^cosrf 

sin"*  ^  (cos*)  m  sin  mcfV  (-l)n  :  Y~    G"     (-l)n 

j  n       nm  ^ 

r       /v      }o<    Tsin(*      &</   '         iKr  ^-    ™  l     J 


) .    cos  m  c^ 


(where  the  subscript^   indicates  differentiation  with  respect  too(  ,   i.e. 

Pn   (cos*V=  ^r  p? (cosot>       X 


Nov;  from   (45)   and   (4.6)   we  h*j.ve 
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(49) 


*£  =  (K^D^r)  , 
EP=  (K2+I^)(^r)  , 
H*  =  (K2t^)(^r)  f 


r        r    p    • 


It  therefore  seems  feasible  (from  comparison  of  (4?)  and  (49))  to  seek 
solutions  of  (49)  in  the  form 


n    run 


(50) 


^l^TZ^v   (-Dn 

l       *   ■  ■■  ^n     run  9 

<fw  =  H oi     a-    (-i)» 
p      * —      n    an  *     '      • 

Substituting  the  quantities  on  the  extreme  right  of  (4S)  into  the  left  sides 
of  (49)  and  them's  of  (50)  into  the  right  sides  of  (49)  yields  a  formal 
identity  which  reduces  to  a  relation  involving  only  the  radial  functions.  Thus 
(5D    TZ  0m  =  (K^Dj)(r  H*n  GnB)  -*  J„  :   (K2-r^)(rJno<n) 

which  is  valid  for  all  four  Quantities.  But  from  the  spherical  Bessel's  equation 
(viz.,  the  equation  of  the  radial  function  resulting  from  the  separation  of  the 
variables  in  the  wave  eauation  in  spherical  coordinates)  we  have 

(52)  (K*+  D*.)(jnr)  =  n(n+l)Jn/r 
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Fe  therefore  obtain  from  comparison  of   (51)   and   (52)   that 

(53)  CXn  =  r/n(n+l) 

Hence  the  Hertz  potentials  may  be  expressed  as 
oo       /„   ,  \       n 


)n 


i*'.i  Hi"-1  %±i|  j_  ^r    is=sUp«  (cos9)(-x: 

i       K    frf  n(n  +  l)      Jn    *~       m   (n+m)!     n 

n— J.  m— o 

X<  AP"   (coso<  ).    cos  m(/ -5—     p"   (cos*)   m  sin  n(J\    =  T~"   1      P'    (-l)n 

/       n  *  '       sin  o<.       n  J        —     "•     nm  j 


V;--  hj„p^,(-i)" 


(54) 


CO 


T  "  =  -  2   7"  in_1     fo*  Xj   j       £"   £    m  (n-)v!   P»   (cos  9)  (-1)1 

Xf  BP™   (cosc*)^   cos  m/+  ^-^     Pj  (cosoOu  sin  m<^   =  JZ  in  O"1)"  , 

y "  =  zz  Jn  pM  (-Dm 

'     p         * "n     nm  • 

Fe  are  now  in  a  position  to  solve  the  problem  of  Fig.  2b  subject  to  the 
second  condition  of  (42):  i.e.,  we  consider  the  original  wave  and  its  image 
simultaneously  incident  on  the  sphere  and  seek  to  determine  the  scattered 
field  EB,  H8,  by  first  deriving  the  corresponding  V  • 

8 

We  have 

v ;  +r ; = c  j,  *;.[  (-Dn  -  (-i)"i  --  *  zz  n  <-u»  j„  p^ 

K  n   m  nm   , 

(55) 

r  >v;  =  eubp^[(-i)n (-!>■]  =2 n n <-Dn Jn p^  , 

-1     ^  n   m  "™ 
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( where  the  notation  ^  ^  indicates  a  sum  over  n  and  m  so  that  the 

n   m 

indices  have  odd  pnrity  while  7       y       indicates  thst  the  indices  have  even 

n   m 
parity).  When  these  rlane  waves  are  incident  on  the  sphere  the  corresponding 

scattered  waves  may  be  expressed  in  the  eetne  form  as  (55)  but  with  j  re- 
placed by  ah,  where  the  h  ,  the  spherical  Hankel  functions  of  the  first 
kind,  are  associated  with  outgoing  spherical  waves.  We  therefore  write 

Mr  =V;+v  ;+^;  =  2  n  3  <-Dn  <v  -;**)  p^    , 

(56) 

The  constants  aR,  an,  are  now  chosen  to  satisfy  the  boundary  condition 

(E„  =  EQ  =  0)    .   From  (45)  and  (46)  it  can  be  seen  that  as  far  as  the  dependence 
"         °     r=a 

on  r  is  concerned  we  may  write  E^,  E '  oc  -  D_  CV*-.)     while 

v.   9    r      r 

Kf*   Fe  0C'V'n»  so  th8t  the  restrictions  on  the  "^  's  are  thatf  Dr(^*'r)  =  V"  -  0^   , 
Imposing  these  on  (56)  (i.e.,  replacing  r  by  a  In  D  (V'r)  and  in^f'n   and  then 
equating  the  coefficients  of  the  angular  functions  to  zero)  yields 

(57)  a^  =  -[KaJn  (Ka)]1  /(&%  (Ka)  V      ,  aj  =  -  Jn(ka)/hn(Ka)   , 

where  the  prime  on  the  square  bracket  indicates  differentiation  with  respect  to 

the  argument,  Ka.   (Contrast  these  with  13). 

Consequently  (56)  and  the  field  derivable  by  means  of  (45)  and  (46)  represent 

the  solution  to  the  image  problem  of  Fig.  2b,  and  it  can  be  seen  that  this  is 

also  the  solution  to  the  boss  problem  in  the  region  8  4^/2   satisfying  the 

required  boundary  conditions  at  r  =  a  and  9  -n/2.     The  latter  is  of  course 

obvious  from  (56)  where  it  can  be  seen  that  t"  has  the  same  angular  depen- 

s 

dence  as  ?f.   +V    and  must  therefore   satisfy  the   same  bound-ry  condition 

at  6  =7772  since  this  condition  depends  only  on  ©  or  derivatives  with  respect  to  6. 
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To  see  analytically  that  the  condition  (I^B  *r  ■  0)Q  ^   .   Is  satisfied, 
we  note  froa  (U5)  and  (1*6)  that  (as  far  as  the  dependence  on  0  is  concerned) 
l'  and  %\  oc  V  ■  while  1^  -0,  B^  Oc  DQ  V"  B.  Hence  *r  and  1^  oc  P™J  (cos  0) 

while  lL  ce  4-?  P*^  (cos  ©),  so  that  all  three  quantities  are  sero  at  ©  »^/2 

9 
as  is  required. 

Hen^e  when  a  plane  electromagnetic  *ave  characterized  by  E  as  in  (4-3) 
is  incident  on  the  boss  on  the  infinite  plane  the  scattered  reflected  Hertz 
potentials  are 

(58)  *  r^  (coB*  }«    cos  ■*-  B  "^ —  m  6ln  n<f  J     9 

f  _.  Pj[  (cos*)         ~) 

X4BP^   (cosol  ),    cos  my+  A  m  sin  m  UP  > 

I     n       *  "     sin  <X  ^  S    • 

The  scattered  reflected  electromagnetic  field  is  then 

(59)  Er  =  EP+  Es 

where  by  means  of  equations  (45)  »nd  (46) 

Es  =  E'8-f  E"B 


(59a) 


t   ^r2iw  s        Ttr3r}e   VTs   ^  sine   3^  J    T  lr  sine^r}^8   J      ^    JTi  * 


Making  use  of  (58)  we  obtain 
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Hn(n  +  1)  a^  hn   ^  M  P(t)  [  A  P(a^  C  -  B  i|i|i  n  s] 


Is  =  IfI_ 


(59*) 


«     _ 

-  i  a  h   V~ 

n  n  *— 

D 


■  «»J,f»  «->- ■    ♦  a  ^ .  ,]|  q 


*aere,for  brevity,      H  =  21"11"1  (2n  *  l)/n(n  +  l) 

M  «  €  (n-«).'/(n  ♦  m)« 

P(t)  »  P°  (cos  ©) 
n 

P(a)  ■  P°  (cob<*  ) 
& 

C  "    COS  B<^ 

S  =  sin  n  i/ 


Similarly 


Hr  =  Hp  ■»■  H8 
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(59c) 


-         -         -  dr  B  •*    sin '6    )<f        r  ^  rj8   ^  '  s   'J 

) 

Making  use  of  (58)  again,  ve  obtain 

H8  =  -^H  Mn  +  1)   «;  ^  H  H  P(t)[   B  P(e)<<C  ♦     A^  .     S  ]? 
n  in 

-wTZNfia'   h     ~T     M    ^^r]    B     ^t   a?  C     ♦     A  F(aLn     si 
/  *—-      )      n     n  — *  sin  9L         sin  ot  •*         -I 

n         (.  m 

"fKrhnT     H     «     P(t)[BP(aUC       *       A^m     Si] 


(59<D 


n       Kr 


6 


+*?F"N(ia'    h     "?M     P(t)    Pa     P(a),  C     -     i   M.  i     si 


■  [Krhn]' 

-a — 

n       Kr 


sin  OL 


a     P(e)        2 
sine* 


C      -     BPU^m     S]^     , 
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*here  Ep  end  HP  are  given   by   eouation   (LA)        Th*  „* 

.>        u.^ion   (U).      The  above   constitute  the  solution 

to  the  problem. 

For  Kr»  1  and  »  n,   we  can  Eubstltute  the  asympt^lc  ^  ^  ^  ^^ 
functions  into  the  above  series,    i.e.,   h^     ^  eiKr,  f^j.^  ^     < 
Also     if  Ka<  1  the  series  ,re  repidlv  convergent  and  may  be  approximated  by 
their^irst  fe«  terms  because  of  the  nature  «f  the  dependence  of  the  coefficients 
an,  an,   on  Ka,    i.e.,  a'  =  12(Ka)3/3   ,     a'  =  i  f^lV  ln±AL(Ka)2n+ X 

1  n  l(2n)lh  (2n-$lTZ~~~       ; 

a;  =  -i(Ka)3/3     ,        a;  =  -ir2lnil2     JKa)2^1 

(See  reference  1,  Appendix  IV.) 

Hence  to  the  order  of  (K8,3/Kr  the  scattered  gutties  MJ  be  approxi- 


mated  by 

o .a. » 


f  g  =  -  i2Ka3A  cos  6  sin  o(  eiKr/r 
(60) 


o 
and 


Ts  =  l^Ka^f  B  sin  9  cos  c^  cos* -M  sin  G  sin  J]  eiKr/r      - 

<A3/r)£[A(2  sin  8  sinrf+co.,,  )  -B  sin  yrco.*  J  "J 
-["'  =°s  e  sin<f+8  cos  e  cos/cos.<7£? 


"E6  =  -  elltr 


•her.  the  signifies  of  the  superscrl pt  „  „  as  ]n  ^  ^   ^ 

For  the  pl.n.  of  incidence,  the  „-pl5n<>,  we  M„  sl„pllfy  ^  ^^  fcy 
con.id.rir,  ^  0  and  8  (hith.rto  considered  as  elMy6  ™ltlw)   as  ^^  only 
for  „)0  end  negative  for  x<0.     We  then  obtain 
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;  ^n    ^oZ    -iKr  cos  (e  '*}    j 
Ei  =[-A[8ln  (9  -«)^  +  cos   (9  -*>   9>B*S  e 

U3')~  ^2    _1Kr  cob  (e  -*). 

H1  =7^B[sir    (9  .^)^cob   (6  -<0   0>    AC^f      e 


iKr  cob    (9+ot  ) 
e  1 


EP  =   £[aln   (6  -NO  ^+ cos   (9+-  )  ^     -*?3 

(U?)  /vn       -?    iKr  cos  (e4<<)  ; 

BP  =  ^B[sin  (9  +  *  r?«  cos   (9  -rrf  )  -8]  +  A^  $   • 

(K2  a3/r)U  (1+2  sin*  sin  9)  S    -B  cos*  cos  9  (<  J  , 


iKr 
0£B  =  _  e 


0n  adding  the  p  and   s   tarns  *•   obtain 

-  .     iKr  cos   (9+-OC  A  sin  {QW  )^A[Cos(9  +  ^)-Qe  e1^  -B[l-Q<    ^1%}  y 
tr   —  e  j 

(6l"r  -^    iKr  COB   (0+/M  B  Bin(94*y?     b[cos(9+*  )   -O^e1^^  A^  ^J  Q  \    \ 

where 

q  =     2Kr  sin2  [  (9  +  *)  /2]        , 

0     =   (K2  aVr)((l+2  Bin*. In  eM^ae)      ...J    , 
q     =    (K2  a3/r)  {  «>e<<.  cos  et  C^'tKa)2    ...        ^    • 

(The  aotationStea)2,  OfW.  -rely  indicates  t*t  the  next  order  ter.s 
ln  the  curly  breeds  9re  to  toe  order  of   (Ke)2  or  that  the  next  order  ter,s 

in  C  are   to  the  order  of   (Ka)    ). 

Tnese  expressions  can  he  readily  roared  1th  those  obtained   for  the 
cylindrical  cose  for  incidence  perpendiculer  to  the  axis  of  the  boss 
(5S6M7M7-).     This  comperleon,  as  «11  as  further  considerations  of  the 
implicative  of     he  abo,e,  will  he  reserved  for  the  lest  section  of  this 
report. 
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B.     Distributions  of  Bosses  . 

Subject  to  essential!,  the  ™  restrictions  as  for  the  cylindrical 
case  the  asymptotic  aolutlon  foe  the  single  W  will  „„.  be  .xtended  t0  obtei„ 
approve  solutions  to  some  problem,  of  scattered  reflection  from  certain 
planar  distributions  of  such  bosses.     The  interactions  of  the  bosses  mil!  be 
neglected  as  previously  and  the  finite  distributions  considered  will  be  such 
that  their  „,M  extent  is  much  less  then  r  sin  9   (or  less  than  r  if  «.  .  0) . 

The  procedure  to  be  followed  In  this  aection  dirf.rB  fro„  that  ^^ 
for  the  -nalogous  cylindrical  problem  in  that  fro.  this  point  on  attention 
•ill  be  reacted  to  the  field  vectors  and  ..  .U1  no  longer  be  concerned 
with  the  Bert,  potentials.    (The  reason  for  this   is  that  it  is  the  simpler 
procedure  to  follow.     The  potential  for  a  distribution  of  scicyllnders  was 
a  si»ple  vector  puantity  having  the  same  direction  at  a   field  point  as  did   the 

"""^  P°tentlElB  "  «*■"*»*.   1«  th.  unit  vector S  .as  comnon  to 
all  the  euantities.     The  potential  for  a  distribution  of  hemispheres, 
consisting  of  a  superposition  of  elementary  potentials  each  charactered  by 
a  different  unit  radial  vector  to  a   field  point-wouK,  necessitate  the  manipu- 
lation  of  entities  possessing  counts  along  all  three   coordinate  vectors  and 
would  result  in  no   simplification.) 

1)     Patterned  Dietributions  . 

Consider  J«    incident  on  5Ntl  -*  bosses  of  radius-a  whose  centers  „re 
spsced  symmetrically  (with  spacing-b)   along  a  line  in  the  xy.nlane  through  the 
origin  making  an  enRleH„Ith  the  x^u   (i  „rjy„  ^  tesseE)_     ^  sMtt^ 

components  of  the   field  reflected  from  the  ray  can  he  expressed  as 

N 

(62)  pS  a  iZZ         qFs  e~lqb  cos(*sin<< 

0-  -N  * 
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,here  the  exponent!.!  ^se  factor  is  chosen  so  thet  the  phase  ia  sero  at  the 
orifln  of  the  cent,,,  toss  <the  reference  origin)  o  -   0,  sea  ehere  sublet  to 

^     *      +**   cp8  .F  =  E.  H)  are  obtained  from  (61)  by 
the  single  scattering  hypothesis  the  F  (F   &,  «; 

,   4      fl^r   9  ,</„,  the  coordinates  with  reference  to  the  q'th 
replacing  r,  9,  cf,  by  rq,  °q>^q> 

*„~e.   thP  procedure  we  will  follow  is 
boss  and  similarly  for  the  unit  vectors.  The  rrocedur 

i   ttt  Rl  except  thPt  for  simplicity  we  rill 
essentially  that  of  reference  1,  III  Bl,  except- 

«  <„  (£,"[}    rather  than  with  the  peneral 
work  with  the  asymptotic  solution  as  in  (61)  rather 

solution. 

Now  if^b«r  sin  9     we  have  approximately 

r  oc  r  -  c :b  sin  8  coe   (^/-^)        , 

q 
(63)  cos  9  -  cos  8q«  -  3>_  sin  9  cos  (f  -<)   cos  9Q  ->  6a*9      , 

<<,*  /+    a^9fi  +  «  ■ 

Haying  the  Utter  approximations  into  the  V    fleeting  ph  sin  9  cos   <,-*>«. 
ln  the  denominator  =nd  similarly  aoproximating  Sq  by  IS  and*,  by«  )  ana  evalu- 
ating the  resulting  geometrical  progression    yields,   essential,  es  for  the 
cylindrical  esse,   tlmt  f  =  °I6  sin^r/slnP     and 
(64)         f  =  FP+  "Is  sin  'W/'/sin  P   ,      P  =  Kb[sin  *  cos  f,  t  sin  6  cos  ( ft  -  <<  fl 1*       ■ 

Por  any  distribution  of  tosses  in  the  xy-pl.ne  .hich  can  V*  decomposed 
into  lines  at  arbitrary  angles  f^rtth  the  x-.xis,  each  line  consisting  of  s 
superposition  of  J  rays  each  charactered  hy  a  different  number  of  tosses   <*  >  . 
w  .  different  spacing   (»).   snd  different  toss  redius   (a),  -nd  al!  lines  and 
rays  having  the  centrsl  element  in  common,  ee  eould  oroceed  by  summing  over 
j  andU   and  obtain  (63)   but  containing 
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(See  discussion  of  expression  21). 

When  I     is  incident  on  a  reflection  lattice  of  bosses  symmetrically  dis- 
tributed with  respect  to  perpendicular  axes  x',  y',    (which  hare  been  rotated  through 
an  angle  ^    from  the  reference  xt  y  axes)  characterized  by??1  b',  Vft  b',  we  obtain  in 
analogy  to  (62) 

(65)  f*  =  tl     r~_    qx,qyF8  e-iKecos(y-+a )  einc< 

whereO  =  J"(qzf  b^,)  ■+  (qy,  by,)2J*  and  tanp'  =  o       b  f/qx,  b^,  ,  and  where 
the  approximations  to  be  employed  are  of  the  same  form  as  in  (63)  but  with  k*  re- 
placed by  jui  '48  and  qb  by  p.  Hence  oroceeding  as  in  reference  1,111  b  7, 
yield 8 

fr  =  FP+°F«  sin/tfx,pxt    Sin^7y,ryl/sinPx,    sinPy,      % 

(66)  n 

I    xi   =  *Vl  coe^sinj-f    sin  0  cos(^-g)J  /? 

P  yi   =  Kby,[  -  sin0sino<-t  sin  6  sin   (^-^U  /2 

For  a  planar  superposition  of  such  distributions  with  a  comron  central 
boss  we  would   obtain   (64.)   but  containing 

y-    r-    3      r  8ln^M(?)Px.3(e)sl"^y3(e)Py.3(g)  -1 
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2)      Uniform  Random  Distributions  . 

For  E     incident  on  a  random  distribution  of  bosses   in   the  xy-plane 
characterized  by  a  density  function  K  which  is   independent  of  the   coordinates 
so   th^t  W  dx  dy  is   the  number  of  bosses   in  ares   dx  dy  we  proceed   in  analogy 
to  the  above  by  evaluating 

(67)  £■  z  mJTPpb  e-ttpeo8^«in*    dx  dy 

where  p  =  (x2^  y2)?  ,   tan/4  =  y/x  . 

(a)   Small  Finite  Distributions 

For  finite  distributions  of  maximum  extent  «r  sin  9  we  employ  essentially 
the  same  approximations  as  in  (63)  and  the  resulting  simple  integrals  are 
readily  evaluated  as  in  reference  1,  III  C  1. 

For  a  ray  of  length  d  by  unit  width  we  obtain 

(6*)  f  .  fP+OjB     Md  SinP/n 

(Pas  in  (63)  with  b  replaced  by  d)   . 

For  a  rectangular  patch  dxf  by  dy, ,  rotated  an  angle  £  from  the  x,  y,  axes 
as  for  the  lattice,  we  obtain 

(69)  Fr  =  £P+  °FS  Md  ,  d   sin  P   sin  P   /  n   H 

x*      y»  x'  '     y«     I     x'  '     y» 

('   x»»Pvi    »   as  in   (66)   with  the  b's  replaced  by  d's) 

For  a  circular  patch  of  radius  d     we  obtain   for  the   plane  of  incidence 

(TO)         f  =  FP+°FE     M7Td2[2  ^    (  ro)/rol  »    «V  K  ^(sir^-f  sin  «)    § 
where  9  is  ±  for  x  £  0. 
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These  results  may  also  be  extended  by  summation  to  a  coplanar  super- 
position of  such  distributions  each  characterized  by  «  different  density 
function,  extent,  and  boss  radius. 

It  can  be  seen  that  at  the  specular  angle  all  of  the  above  expresrions 
for  the  small  finite  distributions  (6£,  66,  68,  69,  "TO) — as  well  as  those 
for  the  superposed  distributions,  provided  that  all  a,  -   a — can  be  expressed  as 


lj 


(71) 


(f   =  fP+  Of8<9?  ) 


0  =g*,<^  =  7T 


V 


where^?  is  the  total  number  of  bosses  in  the  distribution.  Fe  may  therefore 
conclude  that  this  is  also  true,  subject  to  the  initial  restrictions  and  to  the 
same  order  of  approximation,  for  any  distribution  (of  maximum  extent  «  r  sin  9) 
of  bosses  of  equal  radii.   (,See  discussion  of  expression  22). 

(b)  Infinite  Distribution. 

For  an  infinite  distribution  in  the  xy  plane  we  evaluate  (67)  by 
the  method  of  stationary  phase  at  an  arbitrary  field  point,  z  =  D  =  r  cos  9, 
x  =  C  =  r  sin  9  cos  ,  y  ■  y  •  the  phase  factor  in  the  integral  being  replaced 
by  exp  (C  -fp  cosm)  to  correspond  to  the  sketch  below. 
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Purauing  the  alternate  procedure  for  the  cylindrical  ease,  we  write 

the  scalar  coefficients  E,,  E^,  (whose  values  may  be  obtained  from  comparison 
with  (61))  being  functions  of  the  variables  r',  9',^  '  the  coordinates  from  a 
point  on  the  distribution  to  Mie  field  point  P.  We  now  transform  this  quantity 

and  express  it  in  terms  of  the  base  unit  vectors  of  the  reference  frame 

/%  **   ss 

r,  9,t/.  Thus 

es  -  ?f8  •  -?-?  +  Ces  .  'e'e  +  Pes  •£* 

=   (E19'.^E2^'.'r)^(E1^^^E2t?».'9)'9+  (%£' 4)  +     E^'-^W, 

where  the  scalar  products  may  be  obtained  mo?t  readily  by  expressing  both  primed 
and  unpriraed  unit  vectors  in  terns  of  the  reference  i,  j,  k,  components  by 
means  of  the  transformation  equations  given  previously  with  reference  to  (A3). 
Thus 

IT-*?  =  (co«9,cosi<,i+ cos9'sini^ '^  -  sin9k)»(sin9  cosc/i-f  sin  9  sin^j  -*cos9lc) 
=  cos  9f  sin  9  cos  (*/'-*()   -  sin  9'  cos  9   , 
(74)*?''*  =  -  sin  9  sin  (</'-</)     , 

§••6  -  cos  9'  cos  9  cos  (  cf  »-i^) -^  sin  91  sin  9  , 


§'.<?  =  cos  9'  sin  («f  •-</)    , 
<('•?  =  cos  (^'-^) 
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Consenuently 


?E8  =    £  E1t  cob  9»    sin  9  cos   (<{  '-^  )-  sin  9'    cos  9^  -  E2  sin  9  sin(cf  »-^  )?  'r 


(75) 


-V  4E..  £  cos  9'    cos  9  cos   (if  '-«/ )      sin  9'    sin  d~l  -  Ej  cos  9  sin(<f  '-<^  )? 


t{  E^^     cos  9'      sin  (tf  '-/)  +  E     cos   (tf  '-^  )}  Cf 


4 


r-K2a3   (•iKr'/p,)f(6»,^»)     . 


Fe  can  then  express  (67)  as 


(76)  fS  =  -  K2  a3 


foo     /"oo      iKfr'-fc  +cos«)6inoQ   ,        -1      , 
Ml  e  (r«)      f   (9«,<^«)dxdy 

«/-CO    v/-CO 


which  may  be  rewritten  in  terms  of  9'  and  ^f  =  i/  '  -7T  (as  in  reference 
1,  Til  C?)  as 

(77)   Es=  -K2a3M  eiK°  Eln   f2^^2  eiKD(sece»-tan9»sin*eos/Uf) 

Jojo 


X    X  (Q,»  yw+tf"      )   D  tan  9»    sec  9'   d  9'   d/u 

to  yield  on  evaluation   (the  nhase  being  stationary  st   the   specular  angle, 
0»   =  cK    ,tf=0,   orlf'   -IT  ) 


fs  -       v-2  o3 


iKr,T 


(73)  es  -  -  TC  ^  M  f  (9'   =•(    ,C^»   =7T  )(i2fl7K)    sec*    e  -P  -     j 


which  is  valid  except  for  the  vicinity  of^  =^"/2. 
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The  s^ne  procedure } when  applied  to  obtain  4s,  yields  (73)  but  containing 
f'  which  differs  from  the  orevious  f  only  in  that  £.  is  replaced  by  -  ^?E2  and  F^ 
by/?7E  .  The  auantities  _£'  an^  £  can  then  be  obtained  from  (75)  by  replacing 
9'  by  O.   and  c(  •  by  ff  .   Hence  -ve   obtain  the  scattered  reflected  field 

F1"  -  Ep+  Fs 

Es  =  -i3KVsec  e  ~p  "~«Tacos  2*  (co6«(  sin9cos«f-fsin*cos9)-Bcos2o(>  sinQsin^   r 

-*■  fAcos  2  <*(cosd»cos9co8(^-sino(.sin9)-Bcos2oC  cos9sin  Lf~]  9 

-  £  Acos  Pe^eoseC  sint^-r  B  cos^/.  cos  l{~]  <%  r 

=  -i3KVseco<f  cos  26L  eJ+    coe2*     TPL~] 
(79) 

Hr  -   Kp-f-    Hs 

K2  -  -i^3KVsec<*e  -p  -}[bcos2o<.  (cos<*  sin9cos<^4sin  *cos9HAcos  2*  sinQsin^l 

+  [Bcos2o<>  (cos<A.  cos9cos<f-Rin*  sin9)+AcoB  24  cos9sinc^]9 

-  £  Bcos2^  coso<  8in<^  -Acos  2<><  cos  tfj^f 

=  -i3KVseco<  J  cos2*   j£  +  cos  2<*Hjl 

where  V  =  M  ?7Ta^/3   is   the  volumetric  depnrture   from  the   planej    where  Ep  and 
rtp  '3r^    ns   in   (44.)*  and  where   the   subscripts  H  and  J.   indicate   the   components   n 
and  J^  to  the  plane  of  incidence.    (Note  the   relationship  of  the   s  and  p  auanti- 


T 
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ties  of  these  expressions,  from  which  it  can  be  seen  that  were  either  A 
or  B  eaual  to  zero  then  the  scattered  component  (s)  would  be  directly  pro- 
portional to  the  specular  component  (p).) 

For  the  plane  of  incidence,  considering*/  -   0  and  0  as i  for  x  ^  0, 
the  above  simplifies  to 


if  =  Ep  +  E1 


s 


—  » 


iKrcos(9-W), 


a  1RXCOSI  WT«»,  lr  _  -.  y\~  _        .    •) 

EB  =  -i3KVseoo(  e  '^  Acos  2<t[  8in(9-M>  )r>co8(9+*  JeJ-Bcos2^  t<$     • 

(79») 

Hr  -  5P+   H8, 

H8  -  -i/^3KVsec*e  ^  Bcos2,^  sin(9+ot  )r+cos(9+^  ^J+Acos  2</^\ 

where  Ep  and  HP  are  as  in   {UUr) . 

On  combining  the  p  and  s  terms  we  obtain 

E1*  =  (A(l-i3KV  sec<*cos  2<*  )["sin(9+rf  )  'r-J-  cos  (9+<*  )  ^J 

,  ^\->   iKr  cos  (9+rt) 
(79")         -B  (1-  13KV  cos*)if?  e 

Hr  =/fc(B(l-i3KV  coso^)[sin  (9 +*  )  *T-k  cos  (9+^)^9 J 

a  t-t      i«n     .     o_»  \"\7   iKr  cos  (9+»0  • 
■+  A  (1-  13KV  sec  ot  cos  2<*)y>  e  ) 

which  should  be  compared  with  the  analogous  expressions  for  the  cylindrical 
case,  (30"),  which  they  strongly  resemble.  The  discussion  of  the  cylindrical 
case  is  also  fully  applicable  to  the  above. 
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It  should  be  also  noted  that  Es  and  EP  become  formally  Identical  if 

we  substitute  A'  =  i3KVA  sec  <*  cos  2  c*,  and  B*  =  -13KVB  cosc<  Into  the  former 

expression.   It  can  then  be  seen  that  the  "scattered"  Hertz  potentials  for 

the  infinite  distribution  are  formally  identical  with  those  obtained  previously 

for  EP  provided  that  we  replace  A,  B,  inV''  »'2i   as  in  eouation  (54)  with  the 
-  p    P 

above  A1,  B'.   No  doubt  these  Quantities  could  have  been  obtained  directly 
by  operating  on  the  °'V'',  °^"  of  expressions  (60)  but  as  is  apoarent  this 

8      6 

would  hive  been  a  rather  more  complicated  procedure  than  the  technique 
employed . 

The  same  oualitative  results  as  above  (or  as  for  the  cylindrical  case 
for  ft  -7T/2)  may  also  be  obtained  from  physical  considerations  of  a  very 
general  nature.    Consider  the  reflection  of  a  plane  wave,  E%  from  an  in- 
finite plane  surface  covered  with  an  arbitrary  distribution  of  surfsce  irregu- 
larities. If  the  surface  were  perfectly  smooth  the  reflected  wave  would  be 
the  specularly  reflected  plane  wave,  Ep.   The  wave  reflected  from  the  rough 
surface  ^ould  of  course  not  be  plane  but  would  be  pitted  and  corrugated  in 
an  irregular  manner.   If,  however,  the  irregularities  are  assumed  to  be  small 
compared  to  the  incident  wavelength, the  reflected  wave,  having  a  large  specular 
component,  could  be  considered  as  essentially  plane.   We  might  then  express 
the  reflected  wave  as  Ep  eu  or  as  Ep  (l  +  u);  u  being  a  small  complex  quantity 
indicating  th»t  the  reflected  wave,  although  nlane,  has  undergone  an  amplitude 
and  phase  change  on  reflection.   (More  correctly  re  should  emnloy  f£  e11", 
vP  eu-»-  ,  for  the  plane  of  incidence  to  which  the  discussion  will  be  restricted. 
The  following  should  therefore  be  understood  to  anply  to  each  of  these 
components  separately). 


12.  The  above  argument  will  be  found  in  somewhat  greater  detail  in  V.  Twersky, 
On  the  Theory  of  the  Non-Specular  Reflection  of  Sound,  pp.  11-14,  64-66, 
Doctoral  Dissertation,  New  York  University,  1949- 
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This  ounntity,  u,  can  be  but  a  function  of  the  volume  of  the  ir- 
regularities and  their  density  distribution  on  the  plane — or  essentially  the 
volumetric  departure  from  the  plane,  V — the  incident  wavelength,  nn-'  the  angles 
of  incidence  and  observation.  Hence  we  write  u  -  u(^  ,V,ot ,9)  which  from 
dimensional  considerations  yields  u  =  Y  KVg(o(.,9)  where  ¥  is  a  phase  factor 
and  g  is  an  even  function  of  o(.  and  G. 

(The  choice  of  KV,  rather  than  this  quantity  raised  to  some  power 
is  justified  by  there  being  no  a  priori  physical  basis  for  considering 
that  the  scattered  component  can  be  anything  but  a  function  of  the  volume  of 
the  irregularities.  This  consideration  could  of  course  be  obviated  by  assuming 
a  series  in  powers  of  KV  and  retaining  only  the  first  term  since  *re  have 
assumed  that  7.  »  V.) 

Then  by  means  of  essentially  the  same  physical  considerations  that 
justify  the  method  of  stationary  phase  (i.e.,  that  the  rapid  variation  in 
phase  for  regions  not  in  the  vicinity  of  the  snecular  angle  produces  des- 
tructive interference  and  virtually  complete  cancellation  between  the 
waves  from  adjacent  elements)  we  would  surmise  that  practically  all  the 
contributions  to  the  scattered  reflection  "ere  due  to  the  components  at  the 
specular  angle  end  write  u  -^  KVg(  oi)   where  g  is  some  even  function  of  ot  . 
Hence  we  obtain  Er  =  £P(1  -f-yKVg(oO)  (for  either  of  the  components  polarized 
parallel  or  perpendicular  to  the  plane  of  incidence)  which  is  essentially  the 
form  of  the  expressions  (79"),  (30")  derived  previously. 

C.  Distributions  of  Spheres. 

The  results  of  the  preceding  section  can  be  readily  extended  to  the 
problems  of  E1  (as  in  (4-3))  incident  on  similar  distributions  of  whole  spheres. 
The  expressions  for  the  Fs  for  the  small  finite  distributions  (as  in  (6iV)-(71)) 
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are  formally  identical  with  those  that  would  be  obtained  for  the  analogous 
distributions  of  spheres  if  it  is  understood  th*t  °FS  represents  the  asymptotic 
solution  for  the  single  sphere  at  the  center  of  the  distribution  for  Kr'Pl,   Ka<  1, 

The  solution  for  the  single  sphere  may  be  readily  obtained  by  what  is 
essentially  a  simplification  of  the  previous  procedure.  Thus  from  (54-)  we  ob- 
tain the  forms  of  j   ,   and  consider  them  to  be  incident  on  the  sphere.  The 
total  potentials  are  then  of  the  form 


r--r^r,--  ZKJ.+HN,) '„(-«" 


where  the  P   are  as  in  (54).  Imposing  the  boundary  condition  at  r  -   a 
nm 


then  yields 


.  "  EI  an  hn  ?n*t'»  •  ^s  =H  «n  ^  Pnm  (-D 


where  the  a's  are  as  in  (57).   The  field  may  be  obtained  from  the  above  my  means 
of  the  formulas  of  (59, a, c). 

The  rest  follows  as  previously  and  to  the  order  of  (Ka)  far   we  obtain 

(80)   V*g  =  ieiKr(KaVr)[  A(-cos  9  sintf-t  sin  9  cos^cos^)  -B  sin  9  sintf] 

1        =-i/»e   (Ka3/2r)£  B(-cos  9  sin*<-f  sin  9  cos*  coscf  )-f  A  sin  9  sin^J  ] 


and  the  field 
(3D 

■+  [A  sinc((2  coba+ cos  9HB(6in  9  sin«<+  cost^cos  9  coso(-f2  cosi^Jtf  > 


°gS_  e1Kr(K2a3/2r)i^-A(2   sin  9   sinc?-f2   costf  cos  9coscM-cos  i^HBsin^(2cos9+cos«Q]  9 


>Hs  =  ^  (_oEs  ^  +  °E8  t/  ) 

/  *<  9     '  • 
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Had  the  Incident  wave  been  |P   (as  in  (44))   the  above '2^*'   would  be  modi- 

8 

fied  only  to  the  extent  of  replacing  9  by  jf-   9  and  multiplying  °  1^  '  by  -1 

s 

and  the  corresponding  changes  In  the  field  could  be  readily  determined  by 

means  of  (59,  a,  c).  It  can  then  be  seen  that  if  these  two  cases  are  added 

(i.e.,  E1  and  Ep  simultaneously  incident  on  the  sphere,  the  identical  image 

problem  treated  for  the  boss)  that  the  terms  in  A  6in  9  cos4cos<f  and 

-B  sin  9  sin(/  in  the'V'  '  would  then  cancel  leaving  only  -2A  cos  9  sin<< 

and  the  -B  cos  9  sin  ot^  terms  in  the  if      cancel  to  leave  A  sin  9  sin<f-+Bsin  9  cosatcosvf 

as  was  of  course  expected.   (See  expressions  (60)).  To  facilitate  the  extension 

of  the  above  to  the  field,  and  for  reference  purposes,  we  will  also  state  the 

field  obtained  for  £P  incident,  i.e., 

°ES  =  e  r(Ka3/2r)l£'-A(2  sin  9  sin/.  -2  cos  9  cos«(  cost/4  cost/  )4Bsin  t/(cos«<  -2coe9)J'9v 
(61  •) 

-r[A  sin«^(co8  9  -2  eose<H  B(cost^cos  9  cos/ -2  cost/ -sin  9  sine*  )fi  I 

It  should  be  noted  from  this  last  expression  that  if  we  put  <A  -   0,  B  =  0  (to 
correspond  to  a  wave  traveling  in  the  positive  z  direction  and  polarized 
parallel  to  x)  we  would  obtain 


V  =  eiKr 


(F^a3/^)  Accost/ (2  cos  9  -1)  "q+.  sin«^(cos  9  -2)tf| 


as  is  also  obtained  by  the  standard  formulation  of  the  problem  of  the  scattering 
of  a  pla.ne  electromagnetic  wave  by  a  perfectly  conducting  sphere.   (This  1e 
included  here  merely  to  point  out  that  (81)  or  (81»)  could  not  have  been  ob- 
tained directly  from  previous  formulations  of  the  problem  of  the  sphere — as  in 
reference  10 — by  a  simple  transformation  of  the  coordinates  such  as  a  rotation 
around  the  y-axis  but  would  also  reouire  the  rotation  of  the  polarization  in 
the  xy -plane) . 
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For  the   Infinite  distribution  and  en   incident  rave  F   ,   we  obtain  by 
the  method  of  stationary  phase  applied   to    (67)   utilizing  the   field  of   (81), 
that 

Esr  =  -i   KM7Ta3  sec</[    (1  +  2  cos  2c<.)    fJ+  (2+ cos  2*  )  |P_  ]      , 
(82) 

HSI"  =  -i  KM7Ta3  sec*[      (2  +  cob  2<*)      hJ+  (1  +  2  cos  2*)   HPX1  f 

is  the  reflected  wave,  and 

ESt  -  1  KM7Ta3  sec*  [  eJ  +  E^"]     =  i  KM7Ta3   sec*  E1 

(83) 

HSt  =  i  K¥7Ta3  seco<,  H1  f 

obtained  by  replacing  6  by7T-  9  in  the  integral 
before  evaluation .or  more  directly  by  obtaining  the  transmitted  rave  for 
Fp  incident  and  then  replacing  the  p  components  by  the  corresponding  i  compon- 
ents in  the  result — since  from  physical  considerations  both  of  the  transmitted 
waves  (or  both  of  the  reflected  waves)  or  rather  their  components  parallel 
and  perpendiculer  to  the  plane  of  incidence  Trust  undergo  the  same  amplitude 
and  phase  change). 

Had  the  incident  wave  been  £p  we  would  but  interchange  the  i  and  p 
components  in  the  above  expressions.   As  previously,  the  addition  of  these  two 
cases  serves  to  check  the  results  derived  for  the  bosses.  Thus  for  z)0  we  pdd 
(82)  (the  reflected  wave  for  E1  incident)  to  (83)  ?ith  F1  replsced  by  Fp  (the 

transmitted  rave  for  Ep  incident)  and  obtain  -1  -2cos  2<X  •+  1  =  -2  cos  2  o< 

2 
for  the  parallel  component  and  -2  -  cos  2<*+  1  =  -  2  cos  X.  for  the  perpendi- 
cular component. as  was  obtained  previously.   (See  expression  (^)). 

It  should  te  noted  that  4M7TaV3  is  the  volume  density  of  spheres  or  es- 
sentially the  volume  of  materiel  per  soup. re  centimeter  of  distribution  in  the 
xy  -ols.ne . 
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5.      Discussion  and_  Comparisons  . 

In  this   section  further  consideration  rill  be  given   to   the  approximate 
solutions  derived  previously.      For  simplicity  and  to   racilitste  comparison 
the  discussion  will  be  restricted  to  the  plane  of  incidence.    In  addition,   for  the  cy- 
lindrical case  we  will  take  (5   -It/2   (incidence  perpendicular  to  the  axis  of 
the  bosses^while   for  the  spherical   case  we  will  adopt  the  convention  that^  =  0 
and  9  is  positive  or  negative  depending  on  x$0.      It  is  then  seen  that  the  cor- 
responding expressions  for  the  two  cases  resemble   each  other  closely. 

Thus  from  (5'),    (6')  end   (43'),  (4V)   we  may  write  for  both  sets  of 
image  waves, 

E1  ={-  Aj>in(9  -«<)r  +  cos(e  -*)*U   B^e"1  Kr  °0S   (°  "  °°    , 

H1  =/>7{B(sin  (e-^)^+.cos(G-o()eJ^A'aJe-i  Kr  C0S   (Q  "  ^ 
(84) 

BP=/A[sin   <a**)?+coa(Q+«)-Sj     -B  *J  e1  Kr  C°8   (9  *"*  }        I 
jjP=^{    B[sin   (9 -h*)  ^4.  cos   (9**    l^fAlje1  Kr  C°8   (9+<) 

(where  "e  =  k  for  the  cylindrical  case  and  %  -  </  for  the  spherical  case).     For 
the  scattered  components  of  the  reflected  field  we  have  for  the  email  finite 
distributions  from  (17«),    (6l»), 

Ec  -  -^eiKr(7Ti/2Kr)"(K8)2[A(H-2  sinasin  9)6  -B2  cosMcos  9    1c  J 

(85) 

IS  =  -^e       (K2a3/r}[A(l  +  2  sinasin  9)  $  -B  coso<cos  9  tf  J 

H     =  ^[-Fe-9+Fe^J     ; 
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( where  A,  which  yield*  H  at  the  specular  angle  of  reflection,  is  a  "A  -function" 
of  type  and  shape  that  characterizes  the  particular  distribution— see  reference  1, 
p.  29).   The  "scattered  conn>onents"  for  the  infinite  distributions  (henceforth 
to  be  indicated  by  th^  subscript  co)  as  in  (30' )»  C9*)  can  be  written  as 

r  c  r-  ^  a,      «  ai  iKrcos(0  fo() 

E   =  -i  2KVC  sec<*^  A  cos  2<rf[sin(e-foC  ) /r-fcos(9  4»<  )9j-B2  cos^*kj« 


iKrcos(9+«) 

e 

J 


Es     =-i  3KVB  sec*  j    A  cos  2o(£sin(9+«*  ) 'r-rcosO +<*  )0j-B  cos2*  «?J 

(86) 

Hc     =  -i/*2KVc  sec(£  B  2cos2o<  [  sin(9  +<*  )  'r+  cos(9  +*  )0J+A  cos2<*  1c  ^  e" 


3 

iKrcos(9+*) 


Hs     =  -i/fc3KVg  secV^B  cos2W  £    ein(9-r*<   )  "r  +  cos(9  +<*  )eJ-tA  coB2<*/?e 

g 

The  expressions  for  the  reflected  fields,  Fr  =  Fp-f-  F  ,  consist  essentially 

of  two  parte.   The  first  (or  specularly  reflected  plane  wRve  component)  indi- 
cates the  effect  of  the  plane,  while  the  second  (scattered  cylindrical  or  spherical 
waves  for  the  single  boss  and  its  finite  distributions,  or  plane  waves  Tor  the 
infinite  distributions)  indicates  the  effect  of  the  bosses.   It  can  be  seen 
that  as  either  a  or  V  -^  0  or^,  — *  oo  (also  as  r  — ^  oo  for  the  finite  distri- 
butions) then  Fr  — ^  Fp  as  would  be  expected  from  physical  considerations. 

To  the  order  of  approximation  of  the  nbove  exprersions  it  can  be  seen 
that  the  reflected  radial  component  is  completely  specular  while  the  reflected 
9  and  e  components  oossess  both  snecular  end  non-speculer  components  for  the 
finite  distributions  and  only  specular  components  for  the  infinite  distributions. 
For  ell  cases,  however,  the  9  and  e  components  are  reflected  with  altered 
amplitude  and  phase  (as  is  also  the  r  component  for  the  infinite  distributions) 
so  th-:t  the  polarization  of  the  reflected  wave  differs  from  that  of  a  plane 
wave  undergoing  specular  reflection  at  a  smooth  surface. 
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The  angular  dependence  of  the  amplitudes  of  the  scattered  fields  for  the 
analogous  cylindrical  and  spherical  cases  resemble  each  other  closely  as  can 
be  seen  fron  (85)  and  differ  essentially  only  in  the  factor  2  possessed  by  the 
k  tern.  Jor  the  infinite  distributions,  the  overall  resemblance  of  the  cor- 
responding expressions,  although  at  first  startling,  might  conceivably  have 
been  predicted  from  statistical  considerations  of  both  rough  and  striated  sur- 
faces and  the  similarity  of  their  cross-sections  in  the  plane  of  incidence— 
since  a  rough  surface  may  be  approximated  by  just  such  a  distribution  of  hemis- 
pherical bosses  while  a  striated  surface  may  be  approximated  by  the  analogous 
distribution  of  parallel  semi cylinders. 

Tor  brevity  the  reflected  fields  for  the  finite  distributions  may  be 
written  as 

/«  eiKrCO8(et0C){Asin(e-W)^Ai:cos(e^)-0eelqj  -^l-o/°3  ?} 
(87) 
jf.  >\eiKrCOs(9+°(){B8ir1(e^)r+B[cos(9^)-0ee1^+A[l-CQeiqJ  e}   . 

where  from  (17"),  (6l") 

q.  =  2Kr  sin2  [  0+o<  )/2  J  -f  7T  /U 

Q°  -A  (777?Kr)^   (Ka)2|(H-2  sir.  *  sir  9)  +0^   (Ka)2   ...|    } 

Q?  =A  (7T/2Kr)*    (Ka)2j  2   cob  4  cos  ©+^c   (Kn)2      ...  £ 

(88) 

q     =  2Kr  sin2  [r©  +  oO/al     , 


S 


=  A  (K2a3/r)  \  (1+2  sir.  o(  sir.  6)  +  #"8   (Ka)2    ...J 
Qe  -  A  {¥?e?/r)S  cos  <* cos  6+^    (Ks)'    ..,? 
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while  for  the  infinite  distributions 


f=  e1Kr  C08(9^^  A(l-iK1)[8in(e+*^+coe(e+^)^3-BCl-iK2l^ 

(89) 

y.     .      iKrcos(0+oi)c  ^_  -3 

^  <V7e  J  B(l-iK2)  [sin(0^)r4  cos (9+/  )9j-t  A[l-iKj  ^   • 

where   from   (30"),    (7)n) 

K£  «  2  K  Vc   seco([cos  2<-+  C^     (Ka)*    ...J 

K2  =  A  K  Vc   sec*[  c°s2o(  ■*■  (^c    (Ka)2   ...J 
(90) 

K?  =  3  K  V     sec*[  cos  2^-f^     (Ka)2   ...1    . 

1  8  8  * 

K^  z  3  K  V     secotTcos2*    -ftf''    (Ka)2    ...~) 

2  s     L  s         J      . 

The  shove  forms  will  be  employed  to  facilitate  the  following  discussion. 

Ve  will  now  derive  expressions  for  the  electromagnetic  quantities  which 
are  usually  measured  and  for  the  ratios  of  the  components  of  these  quantities 
nolsrized  parallel  and  perpendicular  to  the  plane  of  incidence.  No  detailed 
calculations  will  be  attempted,  the  aim  being  but  to  obtain(for  the  limiting 
forms)  approximate  expressions  whose  behavior  is  readily  apparent  and  which 
are  amenable  for  purposes  of  calculation. 

The  reflected  intensity  (aa  measured  with  a  non-directional  bolometer) 
is  proportional  to  J  =  Re(Er«Er*)  so  that  for  the  finite  distributions  we 
have 
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(91)  J  =  Jn+  Jj   =  |A|  [  1+  QQ  -2  CQ  cos(eta(  )cos  q>|B|2[l-»-  0^  -2  Qe  cob  qj 
while  for  the  infinite  distributions 

(92)  Joo  =  Jnoo+  Jl°°  =  ,A,2(1  *  Kl)  *,*?(l  *  *P     ) 

where  the  appropriate  quantities  may  be  found  in  (88)  and  (90).  It  can  readily 

be  seen  that  Q  vanishes  when  sino(sin  9  =  -  l/2;so  that  the  component  of  the 

reflected  radiation  polsri^ed  parallel  to  the  plane  of  incidence  consists 

only  of  the  specular  contribution  from  the  plane.  At  the  specular  angle, 

6  =  -  o<. ,  this  occurs  when  sin2o<*=  l/2  or  oi  =  IT /l*.     For  the  component 

polarized  perpendicular  to  the  plane  of  incidence  it  is  seen  that 

cos  o<  cos  0  must  eoual  zero  (or  «*.  or  Q  -~ff/2)   for  the  scattered  contribution 

to  vanish.  Similarly  for  the  infinite  distributions  the  scattered  component 

polarized  parallel  to  the  plane  of  incidence  vanishes  f or  <*  -7T  /£.  The 

scattered  component  polarized  perpendicular  to  the  plane  of  incidence  cannot 

legitimately  be  considered  as  vanishing  at  ot  -~7f/2   since  this  region  is 

excluded  from  consideration — see  reference  1,  p.  16,  footnote — but  since  the 

functions,  K2,  do  not  become  infinite  at  0<.  =  7T/2   this  extension  might  be 

accepted  for  A  =  0. 

It  should  be  noted  that  although  the  originally  derived  expressions  were 

accurate  only  to  the  order  of  (Ka)2  0r  (Ka)3  for  the  cylindrical  and  spherical 

cases  respectively  we  can  justify  the  retention  of  both  the  Cr  and  C  terms 

in  the  above  to  only  the  lowest  order  for  moderately  large  vcluer-  of  ^  . 

k2*  ,.  (v   v4 


Thus  we  write  0°  -  ^f  fi(Ka)2-+  f2(Ka)4  ../}    and 


(QC)  =  A  [f^(Ka)4-r2  t^   f2  (Kb)6  ..^ where  f  is  substantially  lesi 


than 
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f     (about  1/2  or  less) —  and  similarly  for  Qs — so  that  we  may  retain  only 
the  terms  to  the  order  of  A  (Ka)2,  A  (Ka)^,    in  Q  and  only  the  terms 
A2(K&)U*  A  2(Ka)    ,    in  Q2  on  the  assumption  that  42f^(Ka)^»  ^f2(Ka)4, 
2^2f?(Ka)     "fr&  f^(Ka)    .      The  same  considerations  may  be  also  applied  to  the 
infinite  distributions  for  M  large. 

The  ratios  of  the  reflected  intensities  polarized  parallel  and   perpendi- 
cular to  the  plane  of  incidence  may  be  expreseed  as 

j                |AlZ[lf  Qq  -2  Qq  cos(9-U)cosc|3 
p     =     _JL  = f 

JJ-  |B|2[     1+  Q^  -  2  Qe  cos  q  ] 

(93) 

0  «,   -_  JL00  =      Ul2(l-»Kf) 
^°°        JJ.oo  |Bi2(H-K^) 

where  the  appropriate  ouantities  are  defined  in   (88),    (90). 

We  will  consider  only  the  limiting  forms  of  these  expressions.     Thus 
for  the  finite  distributions,   if  0«1  we  neglect  02  and  obtain 

(>      =   (|A|2/|B|2)^1-[Q9  cos   (6+<<)    -  0el  2  cos  q^ 
O        =   (  |A|2/|B|2)$  l-.&[(l+2  sin*  sin  6)   cos(6+^  )-2  costfcose] 
X   (2  7T/Kr)^(Ka)2  cos6  Kr  sin2[(o(+e)/2j fTT/^C       f 
P    s     z   (  |Al2/|B)2){  1- ^["(14  2  Bin*  sin  6)   cos(6 +o<)   _  cos  <*  cos  e] 
X     (2K2a3/r)    cos   (2  Kr  sin2[  (ol+  9)/2j)? 
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which  ere  functions  of  all  the  parameters  of  the  problem  and  possess  extreme 
depending  on  Kr,  8,  and  OC.   (To  plot  these  quantities  as  functions  of  9  it 
should  be  noted  that  the  phase  term  is  verv  rapidly  oscillating,  since  KrV)  1, 
so  th8t  it  would  suffice  to  consider  only  its  limiting  values,  ±:  1.  Similarly 
the^i  -function,  except  for  its  first  maximum  at  9  =  -  o^,  may  be  neglected 
for  an  approximate  plot).  To  this  order  of  approximation  it  can  be  seen 
that  in  addition  to  the  vanishing  of  the  scattered  components  as  determined 
by  the  C's  mentioned  previously,  the  scattered  contribution  parallel  to 
the  plane  of  incidence  also  vanishes  at  9  -fo(  =  7T/2. 

At  the  specular  angle,  9  =  -(A  and  A  — >  /??     so  that  the  above  reduce 
to 


(95) 


pc  =  (|A|2/|bI2)[  l-rWOTAr)*  (*a)2J  , 
P  8  =  ()A(2/|B|2)[l-r^sin2c*  2  K2  a3/rj  ^ 


Both  of  these  ire  explicit  functions  of  the  number  and  size  of  the  bosses, 
the  wavelength,  and  the  distance  of  observation  .  Only  the  latter,  however,  de- 
pends on  the  angle  of  incidence  and  has  i is  maximum  value  at  <*.z.  Tf  /2  and  its 
minimum  value,  which  is  also  its  specular  value,  at  <*  -  0.  This  is  a  signi- 
ficant difference  between  the  behavior  for  the  cylindrical  and  spherical 
cases  and  its  experimental  verification  should  be  attempted. 

These  values  also  indicate  that  the  maximum  values  of  the  O  *s  as 
functions  of  9  occur  at  essentially  the  specular  angle|  primarily  because  of 
the  behavior  of  the  A  -function. 

For  the  case  Q2  >^  1,  as  for  large  A  but  with  Ka<<;  1,  we  obtain,  on  the 
assumption  that  the  remaining  terms  may  be  neglected,  that 
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2,,_         _x2 


p   =  (lAl2/|Bl2)  of/Q2 
(96)     n       -   (|A|2/|B|2)  (It  2  sin*  sin  6)*/(2  cos*  cos  9) 
=  (|AJ2/lB|2)  (1+  2  sin«Csin  9)2/(cos  xcos  9)* 

fit  should  be  noted  that  the  negative  sauare  root  of  the  above  O         for  an 
unpolarized  incident  wave  (A  -  B  -   I),vit.,-(l-r  2  sinotsin  0/2  cos oi  cos  9), 
was  the  only  quantity  considered  by  Rayleigh  with  the  assumption  that  the 
specularly  reflected  components  (due  to  E  )  could  be  screened  off — an 
assumption  that  might  be  justified  except  in  the  vicinity  of  the  specular  angle. 
Rayleigh  found  oualitative  agreement  between  the  above  and  experimental  re- 
sults obtained  for  ruled  metallic  surfaces  or  surfaces  covered  with  raised 
linear  protuberances  for  the  cases  <XvO,  Q~Tf/2   or  9~0,  e<  -~W/2   and  for 
the  case  where  oi  and  9  both  made  abtuse  angles  on  the  same  side  of  the 
surface  normal.  He  did  not,  however,  find  experimentally  that  the  parallel 
component  of  the  reflected  radiation  vanishes  as  1+2  sin  o^sin  9  approached 
zero  and  surmised  that  this  might  be  due  to  the  fact  that  the  radii  of  the  pro- 
tuberances were  not  small  compared  to  the  wavelength,  i.e.,  that  the  condition 
Ka<l  had  not  been  satisfied.  This  could,  of  course,  also  be  due  to  the 
fact  th8t  the  neglected  specular  contribution  was  significant  for  the  angles 
investigated.}  It  can  be  seen  that  these  auantities  (96)  depend  essentially 
only  on  the  nngles  of  incidence  and  observation  and  that  the  ration  of Pc  to 
P  fl  is  lA;  or  the  effect  of  the  hemispherical  boss  in  introducing  a  relative 
change  of  polarization  is  four  tim^e  as  great  as  for  the  semi cylindrical  boss. 
It  car.  also  be  seen  that  the  reflected  radiation  is  polarized  perpendicular 
to  the  plane  of  incidence  when  1+2  sinXein  9  =  0  and  polarized  parallel  to 
the  plane  of  incidence  when  cos^cos  9-0.  At  the  specular  angle  the  above 
reduce  to 
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(>c  =   (|a|2/|b|2)   cos2  2d  A  cosV      , 

(97) 

p   =  (|A|2/|B|2)  cob2  Ztf/cos1*  <X    ' 

so  that  the  parallel  component  vanishes  when  cos  2oC  =  0  or  A  =  TC/U  and  the 
perpendicular  component  vanishes  when  eos(X=  Ooro^:  7T/2.   (In  view  of  our 
comments  on  fiayleigh's  procedure  it  should  be  remembered  that  the  specular 
contributions  have  been  ignored  in  the  statement  of  (96),  (97)  so  that  when  the 
scattered  components  become  small  the  complete  expressions  as  in  (93)  or 
those  of  (94)  should  be  employed.  Consequently  the  vanishing  of  a  scattered 
component  still  leaves  the  specular  component) . 

For  the  infinite  distributions  the  second  expression  of  (93)  gives  some 
indication  of  the  orders  of  magnitude  to  be  expected.  It  can  be  seen  that 
for  K2»l,  as  for  large  M  but  with  Ka<<  1,  the  expression  yields  those  of  (97) 
which  are  valid  only  subject  to  the  above  qualifications. 

The  real  part  of  the  radial  component  of  the  complex  Poynting  vector, 
the  average  outward  energy  flow  as  measured  with  n   directional  device,  can 
be  expressed  as 

(93)       S  «  fie  (ij  Bj*  -  ^  Ip/2  -  S,  +  Sx      # 

Hence  for  the  finite  distributions 

S  =   {y/2)  £  |A|2    (cos   (9+wHQ2  -  Qq  cos  qf  ltcos   (9+rff]) 
+  |Bl2    (cos   (9+rf  )+  C^  -  0     cos  q[  1-rcos    (9+o<   )    ")? 

while   for  the  infinite  distributions 
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(100)         S^   =  (y/2)    cos  (6+rf  )£  |A|2(1+K^)+-|B»2(1+K^)j 

where  the  quantities  are  as  defined  in  (38),  (90).  The  behavior  of  these  ex- 
pressions as  regards  their  dependence  on  the  Q's  and  K*s  is  the  same  as  that 
of  the  J's  of  (91),  (92)  and  in  addition  the  specular  components  vanish 
when  (9+d  )   =  T/2. 

The  ratios  of  the  components  polarized  parallel  and  perpendicular  to  the 
plane  of  incidence  can  be  expressed  as 

S     |a|2(cos(9  -»•(*  )+  Qq  -  Qg  cos  q  [  1+  cos(9  r <*   )]  ) 


3J-    |B|2(co8(9  +  e(.)+Q^  -  Q  cos  q[l+  cos(9  +  <*  )Q) 


for  the  finite  distributions,  while  for  the  infinite  distributions  we  obtain 

^oo  =  ?oo  a*  *«  (93). 

For  the  case  Q  <if  1  we  obtain 

IT    =   (|Al2/|B|2)^  l-(Qft  -  Qe)[  1+sec   (04-oOj    cos  q^    j 

<TC  -  (  |A|2/|B|2)  {  1-A["  1-2  cos(9-ro<  )]  [   W  sec   (9  +otj] 

y    (  ^/2Kr)3(Xa)2  cos   (2  Kr  sin2[  (9-K*)/2]    -r7T/4)$    f 

0%  =  (|Aj2/lB|2)  j  1-AC  1-  cos(e+o<)H  sin*sin  9^1 

X  [1  +sec   (9  -f  *  )1  (K2  a3/r)    cos   (2  Kr  sin2  [  (9  +<*  )/2]  H 

To  this  order  of  approximation  it  can  be  seen  th?t  both  these  quantities  re- 
duce to  their  specular  components  when  1-J-  sec   (9-fo^)   -  0  or  9-ro<  =  #" 
and   that  in  addition  (Tc  also  hfc  s  its  specular  value'  at  1-2  cos   (9  ♦  <<)  ■  0 
or  9+ot  =7T/3.     Now  if   (9  +o(  )  =  7f  /2  these  expressions  reduce  to 
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(T-   (|Al2/(B|2)  QQ/Qe  where  [1  +  2  sinoUin  (7T/2  -otQ  /?  cos*  cos  (7f/2  -  c<) 
«  (1+2  sin<<  cosci  )/2  cosoCsin*  =  cosec  2<rf-tl.   Hence 
rc  -   (|Al2//Bl2)  (cosec  2<*+  1)  -  <Hs/2,  so  that  when  <^-7TU%(J9=  If /U) , 
<7~c  reduces  to  its  specular  value  while  ff"e  is  twice  its  specular  value. 
This  is  another  significant  difference  between  the  cylindrical  and  spherical 
cases  and  its  experimental  verification  should  be  attempted. 

It  should  be  noted  that  the  above  <r's,  for  posit? ve  values  of  the  C'3, 
exhibit  essentially  the  same  behavior  as  the  /o  'a  in  (96)  and  might  therefora  account 
for  the  oualitative  agreement  with  experiment  vifhich  Rayieigh  found  for  this  range 
of  values — the  vanishing  of  the  specular  components  obviating  any  additional 
assumptions  as  to  their  significance.  These  results  are  also  preferable  for 
comparison  with  Rayleigh's  experiments  since  he  employed  directional  visual 
devices. 

At  the  specular  angle  the  above  quantities  (102)  become  identical  with 
those  of  (95). 

For  the  case  CT  »  1  if  we  neglect  the  remaining  Quantities  ps  previously 
the  expressions  for  <T  become  identical  with  those  of  (96)  f or /O  subject  to 
the  same  Qualifications. 

We  will  now  consider  the  total  reflected  intensity  and  energy  flux 
for  the  case  where  the  incident  wave  is  unpolarized.   (The  ratios  of  the  re- 
flected components  polarized  parallel  and  perpendicular  to  the  plane  of  inci- 
dence can  be  readily  obtained  by  setting  A  =  B  =  1  in  the  previous  expressions). 

For  the  finite  distributions,  letting  A  =  B  =  1  in  (91),  (99),  we  obtain 

J  =  2-f  (QQ+-Q2,)  -  2[QQ  cos  (9-M  )+Qe3  cos  q 

(103) 

S  =   (/??/2)   £2  cos   (9+o<  )+    (c|+Q2)    -   (Qe+Qe)[n.cos(8-J-^)    cos  q]  2 
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(104) 


For  the  case  C<<(  1  *e  neglect  the  terms  in  0^  and  obtain 

Jc  =  2  $  1-^f  (14  2  sin*8ine)co6(9+c(  )  +2  coed  cos  9] 

X (  T  /2Kr)2    (Ka)2  cos(2  Kr  sin2  [  (9  -+<x)/2  1  4  7fA)  ( 

J     =  2  ^  1- A[  (142  sino(sin  9)cos(9  -+p<)4  cos  <*cos  9 

X  (K2a3/r)cos(2  Kr  sin2[  (6  •+  oC  ) /2}  )?  , 

S     = /*?{cos(9  +  <X)-<A["    1  +  2  cos(9  -*)J£  l+cos(9-ro03 

X   (7T/8Kr)^   (Ka)2  cos(2  Kr  sin2  [  (9  ■+•<  )/2]  ■+  7T  A^V      , 

bs  =  /??lc°B(Q+<<  )"  At  l+eos(0  -<*)  +  sin*  sin  ejf  1-r  cos  (9 ->■  *<  )J 

X    (K2a3/2r)cos(2  Kr  sin2[  (9  +o()/2j)2 

It  should  be  noted  that  the   litter  expressions  reduce  to  S  =  -   {1h/2)(QAQ e)cos  q 

when  9-f<*  =  7T  /2.      Hence  the  above  simplify  to  Sc  =  -  <ff&  (1+2  Bin  2*0   X 

(  7T  /8Kr)^(Ka)2cos(Kr  +  7rA)   and  S„  =  -,#A(2  +  3  sin  2<* )  (K2a3/4r) cos  Kr  for 

s  7 

9  =7T/2  -tk  . 

At  the  specular  angle  the  expressions  of  (104)  reduce  to 

Sc  =  (y/2)   Jc  =^[1-^(1^2  cos  2<X)(ir  A  Kr)*(Ka)2J 

(105)  - 

Sb  =  (J/2)   J     -nrjl   1-99(1^3  cos  2*)(K2a3/2rO 

These  are  essentially  the  extrema  of  (104)  with  respect  to  9.  It  should 
be  noted  that  wnen  14  2  cos  2e<  =  0  -4  <X.  -  60°  and  when  143  cos  2oC  =  0  —:>  0*^55°, 
the  quantities  reduce  to  their  specular  values — i.e.,  for  these  "critical  angles" 
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the  reflection  at  the  specular  angle  is  completely  specular.     Hence  for  all 
angles  of  Incidence  less   than  60  ,  550,   for  the  cylindrical  and  spherical 
cases  respectively  these  indicate  a  minimum  in  the  intensity  and  energy  flux  at 
the  specular  angle,  while  for  values  greater  than  these  critical  angles,   these 
quantities  have  their  maximal*  at  the  specular  angle.     The  extrema  of  (IO5)  with 
respect  to  <*  occur  at  sin  20*.*  0  or  o(  =  0,    7T/2t—the  former  being  the  value  at 
the  minimum  and  the  latter  at  the  maximum. 

Had  the  incident  wave  been  polarised  parallel  to  the  plane  of  incidence, 
as  obtained  by  setting  A  ■  1,  B  «  0,  in  the  initial  expressions,   then  only  the 
Q_  term  would  be  present   in  (lOh)  so  that  the  scattered  component   of  (105)  would 
be  proportional  to  cos  2  c*  and  the  critical  angle  for  all  the  quantities ,  at 
which  the  reflection  at  the  specular  angle  was  completely  specular,  would  occur 
at  cos  2<X»0ora"7T  /U,  the  expressions  having  their  minimum  values  below 
this  value  and  their  maximum  values  above  it.     The  extrema  of  (105)  also  occur 
at  sin  2o(»  0  or  o<  «  0,  T/2  as  above. 

lad  the  incident  wave  been  polarized  perpendicular  to  the  plane  of  incidence, 
as  obtained  by  setting  A  »  0,  B  «  1,   in  the  initial  expressions,   then  only  the  0. 

© 

term  would  be  present  in  (lOU)  so  that  the  scattered  component  of   (105)  would  be 

2 
proportional  to  cos   o<.   and  the  critical  angle  would  be  0/  »7T/2  (i.e.,   the  re- 
flected wave  at  grazing  indicence  would  consist  of  only  the  specular  component)  and 
the  reflected  radiation  would  be  a  minimum  at  the  specular  angle  for  all  other 
angles  of  incidence.     The  extrema  of  (105)  occur  when  cosot     sintf*  0  or <*sO>  77%? 

as  previously. 

2 
For  the  case  Q  »1,  we  have  on  neglecting  the  Q's  that 

J     =  2  +A2[  (1+2  sinotsin  6)2+  (2  cos*  cos  6)' J  IT  K?b/*/2t      , 

3-2  +A2\  (142   sino<sin  9)2-f  (cos^cos  e)'J  K^a6/r2        , 

*/jtf£eos   (9-f  «<)-r  A2[  (l-r2  sin<*sin  9)2-f(2  cos^cos  6)2J  7[  K3aV4rJ  $ 

3g  =  <J?(cos   (9+<*  )+A  2[  (1+2  sin*  sin  9)^+  (cos  o< cos  Q)2~]  K^afyar2  J     ; 


E 

(106) 


&C 
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which,   at  the  specular  angle,   reduce  to 


S     =   (7/2)  Jc  -ip{  1  -f<^2(?  cob2  2<*-t-2  cos  2cX-f  DTTKV/Ar  j    ( 
Ss  =   {/»/2)   Jg  =^?{   l+>?72(5  cos2  2^-f-  2  cos  2<*  41)   K*66/8r2J    , 


These  are  essentially  the  extreme   of   (106)   with  respect  to  0  and  indi- 
cate 8  maximum  at  the  specular  angle  for  all  angles  of  incidence — the  expressions 
not  reducing  to  their  specular  values  for  any  real  values  of  OC .     The  extreme 
of  (107)   occur  when  sin  2  c<  (2  cos  2*4    1)   =  0  or*  =  QtTT/2,7T/3  for  the 
cylindrical  case  end  when  sin  2  o<  (5  cos  2  di  I  1)   :  0  or 

C(~  0,7f /2,  or  t)(^  51°  for  the  spherics]  case — the  expressions  possessing 
maxima  at  the  first  two  values  ofoC  (that  at  crt.  =  0  being  the  greater)  and  a 
minimum,  whose  value  is  still  greater  than  the  specular  value,  at  the  third 
value  of  o{ . 

Had  the  incidsat  vave  been  polarized  parallel  to  the  plane  of  incidence, 
then   (107)   would  contain  cos     2o<.  so  that  the  expressions  would  reduce  to 
their  specular  values  ato<  =  7T  /U,    (as  above  for   (105))   but  would  be  the  maxima 
of   (106)   for  all    other  values  of  c<  .     The  extrema  of   (107)   with  respect  to  *( would 
then  occur  et     2  cos  2  <<  sin  2oC  =  0  or  c^~  0,  7T/2,  TT/U',   being  maxima  for  the 
first  two  values  of  o^and  a  minimum   (the  specular  value)   for  the  third. 

Had  the  incident  wave  been  polarized  perpendicular  to  the  plane  of  inci- 
dence then   (107)   would  contain  cos4-0<  so  that  the  expressions  would  reduce 
to  their  specular  value  at<x!  = *fl/2   (as  above  for  (105))   but  would  be  the 
maxima   for  ell  other  values  of  oi.  .      The  extrema  of   (107)    then  occur  when 
4  cos^o<    sinX  =  0  orc<  =  0,  7T/2,   being  a  maximum  for  the  first  value  of  oi 
and  a  minimum  (the  specular  value)    for  the  second  value. 

For  the  infinite  distributions  we  obtain  on  setting  A  =  B  =  3    in 
(92),    (100)   that 
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J       -2  +  - 
oo 


(10ft) 


Hence 


Sco  =  ^^   cos(9+*)(2  4K^+  K2^)  . 


Sc  =  (^?/2)cos(e-4-^)Jc  =^coa(e-K)  |l-*  2  K2^  sec2*  (2  cos2?* +2  cos  2<*-+l)J  , 
(109) 

S  =  (^  /2)cos(e+o()Je  =^cos(9  +  <<)  [jL+ (9/8)!^  sec2*  (5  co622o< -f2cos2  -f  1)]. 

At  the  Epecular  angle,  cos  (6+^)  would  be  replaced  by  1,  to  yield  the  maximum 
values  for  the  S's. 

It  is  unfortumate  that  no  Quantitative  experimental  data  are  at  present 

it- 
available  with  which  the  above  could  be  compered.   The  expressions  indicating  a 

minimum  in  the  radiation  reflected  at  the  specular  angle  until  some  critical 

angle  is  reached  and  a  maximum  in  the  reflected  radiation  above  this  critical 

angle  are  of  particular  interest  end  their  experimental  verification  6hould  be 

attempted.  As  was  mentioned  in  the  introduction,  minima  at  the  speculer  angle 

3 

had  been  observed  by  Skudrzyk  when  sound  underwent  non-specular  reflection  at 

a  rough  surface  (in  the  shape  of  a  finite  rectangle)  whore  irregularities  were 
of  the  same  order  of  magnitude  as  the  wavelength — somewhat  substantiating  the 
expressions  derived  for  the  acoustical  case  which  indicated  essentially  the 
same  behavior  in  this  respect. 

It  should  be  emphasized  that  the  solutions  obtained  for  the  various  distri- 
butions can  only  be  considered  as  first  order  of  scattering,  or  single  scattering 
approximations  in  the  sense  that  only  the  primary  or  plane  wave  excitation  for 
each  boss  has  been  taken  into  account.  Attempts  h?ve  been  made  to  determine 

*  The  writer  has  just  become  aware  of  some  work  being  done  by  A.  Aden  and 
J.  Sivak,  Cruft  Laboratory,  Harvard,  which  might  have  some  bearing  on  the 
above. 
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the  limits  of  validity  of  this  approximation  airi  to  obtain  more  accurate 
solutions.  Thus,  the  problem  of  the  finite  greting  of  parallel  cylinders  has 
been  formulated  to  take  into  account  the  excitation  of  each  boss  caused  by 

the  raves  scattered  (as  a  result  of  the  primary  excitation)  by  its  closest 

13 
neighbors.   It  was  shown  that  when  Kb})  1,  the  largest  correction  terms  are 

l/2 

proportional  to  (Ka)^/(Kb)  '   so  that  they  vanish  more  rapidly  than  the  primary 

2 
excitation  terms  which  are  of  the  order  of  (Ka)  .   (Similarly  it  can  be  inferred 

thut  for  the  spherical  case,  viz.,  a  ray  or  lattice  of  spheres,  that  the 

largest  correction  terms  would  be  proportional  to  (Ka)  /Kb  as  compared  to 

(Ka)   for  the  primary  excitation  terms). 

The  above  formulation,  however,  suffers  from  the  fact  that  the  excitation 

of  only  the  nearest  neighbors  has  been  taken  into  account  when  as  can  be  seen, 

from  the  form  of  the  correction  terms  and  their  slow  convergence  for  moderate 

values  of  Kb,  we  should  consider  the  interactions  of  all  the  elements  of  a 

distribution.  To  remedy  this,  formal  expressions  have  been  derived  for  the 

scattering  of  a  plane  wave  by  any  configuration  of  parallel  cylinders 

which  take  into  account  all  possible  interations  of  the  scattered  radiation — 

i.e.,  a  solution  to  the  scattering  problem  is  obtained  which  simultaneously 

satisfies  the  reouired  boundary  conditions  at  the  surface  of  every  cylinder.   (An 

attempt  is  also  being  made  to  treat  the  analogous  problem  for  a  distribution  of 

spheres).  These  general  expressions  h"ve  been  applied  to  such  problems  as 

the  prating  of  parallel  cylinders  to  obtain  second  order  of  scattering  solutions 

and  will  appear  in  a  subsequent  report. 


13.  V.  Twersky,  J.  Acoust.  £>oc  ^mer.  22  54-0  (1950)  or  see  reference  12 
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APPENDI7 
Extension  to  Absorbent  Bosses  , 

Although  it  has  not  been  found  possible  to  derive  solutions  to  the 
problems  of  the  boss  on  the  infinite  plane  when  both  plane  and  boss  are  composed 
of  absorbent  material,  the  analysis  of  the  present  paper  can  be  extended  to 
consider  the  identical  problems  of  non-specular  reflection  from  absorbent 
bosses  on  a  perfectly  conducting  plane. 

For  such  surfaces  we  would  replace  the  previous  scattering  coefficients 
An,  Ajj,  as  in  (13)  (for  the  case  8  -  "V  /2)   and  a  ,  a  as  in  (57)  with  the 
scattering  coefficients  obtained  for  the  problem  of  the  scattering  of  a  plane 
wave  by  an  absorbent  cylinder  (for  incidence  perpendicular  to  its  axis)  and  by 
an  absorbent  sphere  respectively.  Thus,  where  the  parameters  of  the  absorbent 
obstacle  are  expressed  by  Kc,yL*  and  Kfl»A^8  respectively,  we  obtain  on  im- 
posing the  boundary  conditions  that  the  tangential  components  of  E  and  H  be 
continuous  across  the  surface  of  the  obstacle  that  the  new  scattering  coefficients 
to  be  inserted  in  ecuation  (15*)  and  (59)  respectively  may  be  expressed  as 


A'  "   M:K  Jn(Kca)JA(Ka)-/uKcJA^-ca>Jn(Ka) 
An  "  "  //cK  Jn(Kca)H^(Ka)-//KcJ;(Kca)Hn(Ka) 


yU  CK  J^(Kca) Jn(Ka)-/MKcJn(Kca) Jn(Ka) 
n     //CK  j;(Kca)Hn(Ka)-//KcJn(Kca)Hn(Ka) 


a 

n 


,  .  _  ^2Jn(Ka)I:KgaJn(K6aO'-/,K^n(Ksa)[Ka,1n(Ka)3, 
HK2hn(K8)tKsajn(KEa0,-/W^Jn(KEa)CKahn(Ka)T 

fj  hn(Ka)[Ksajn(KEa)],-^6jn(Ksa)[Kahn(Ka)], 
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It  can  reedily  be  seen  that  as  K  -)  oo^as  is  reouired  for  a  perfect 
conductor,  that  the  A  reduce  to  the  values  of  (13).  Similarly  when  K  -^  oo 

II  6 

(remembering  that  then  Jn(Kga)  0C(K8a)~  )  the  a  reduce  to  the  values  of  (57). 
For  Ka,  Kca<l,  approximate  forms  of  the  above  may  be  obtained  by  means 
of  the  formulas  in  the  appendices  of  reference  1  and  the  asymptotic  expressions 
obtained  for  the  single  boss  can  then  be  readily  extended  to  the  distributions 
di.*- cussed  in  this  paper. 
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